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Abstract. We are concerned with rigorous mathematical analysis of shock diffrac- 
tion by two-dimensional convex cornered wedges in compressible fluid flow gov- 
erned by the nonlinear wave system. This shock diffraction problem can be formu- 
lated as a boundary value problem for second-order nonlinear partial differential 
equations of mixed elliptic-hyperbolic type in an unbounded domain. It can be 
further reformulated as a free boundary problem for nonlinear degenerate elliptic 
equations of second order. We establish a first global theory of existence and 
regularity for this shock diffraction problem. In particular, we establish that the 
optimal regularity for the solution is C°'^ across the degenerate sonic boundary. 
To achieve this, we develop several mathematical ideas and techniques, which are 
also useful for other related problems involving similar analytical difficulties. 



1. Introduction 

We are concerned with rigorous mathematical analysis of shock diffraction by two- 
dimensional cornered wedges whose angles are less than tt in compressible fluid flow 
governed by the nonlinear wave system. The study of shock diffraction problems can 
date back 1950's by the work of Bargman [3], Lighthill [24^ I25j . Fletcher- Weimer- 
Bleakney |14] . and Fletcher- Taub-Bleakney |13| via asymptotic or experimental anal- 
ysis. Also Courant-Friedrichs [lOj and Whitham |27j . 

In this paper, we develop several mathematical ideas and techniques through the 
nonlinear wave system to establish a rigorous theory of existence and regularity of 
solutions to the diffraction problem. The nonlinear wave system consists of three 
conservation laws, which takes the form: 

Pt + + = 0, 
(1.1) mt+Pxi = 0, 

nt + Px2 = 0, 
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for (t,x) E [0,00) X ]R^,x G M^, where p stands for the density, p for the pressure, 
{m,n) for the momenta in the (xi, X2)-coordinates. The pressure-density constitu- 
tive relation is 

(1.2) p{p) = py-f, 7 > 1, 

by scaHng without loss of generality. Then the sonic speed c = c{p) is determined 

by 

Notice that c{p) is a positive, increasing function for all p > 0. 

The two-dimensional nonlinear wave system (jl.ip is derived from the compressible 
isentropic gas dynamics by neglecting the inertial terms, i.e., the quadratic terms in 
the velocity; see Canic-Keyfitz-Kim [5]. Also see Zheng [28] for a related hyperbolic 
system, the pressure gradient system of conservation laws; the same arguments 
developed in this paper can be carried through to establish a corresponding theory 
of existence and regularity for the pressure gradient system. 

Let So be the vertical planar shock in the (t, x)-coordinates, t G M+ := [0, 00), x = 
(xi,X2) £ M^, with the left constant state Ui = (pi,mi,0) and the right state 
Uq = (/9o,0,0), satisfying 

mi = ViPiPi) -PiPo)){pi - Po) > 0, pi> pq. 

When 5*0 passes through a convex cornered wedge: 

W := {{xi,X2) : X2 < 0, —00 < xi < X2 ctan 6^}, 

shock diffraction occurs, where the wedge angle 6^ is between — vr and 0; see Fig. [TJ 
Then the shock diffraction problem can be formulated as the following mathematical 
problem: 




Figure 1. Initial-boundary value problem 



SHOCK DIFFRACTION BY CONVEX CORNERED WEDGES 



3 





Figure 2. Shock diffraction configuration 



{p,m,n)\t=o 



Problem 1 (Initial- Boundary Value Problem). Seek a solution of system 
(jl.ip with the initial condition at t = 0: 
(1.3) 

(pO) 0, 0) in {xi > 0, 2:2 > 0} U {{x2 — xi ta,n9y^)xi > 0,X2 < 0}, 
(pi,mi,0) in {xi < 0, 2:2 > 0}, 

and the slip boundary condition along the wedge boundary dW : 

(1.4) {m,n) -uIqw = 0, 

where v is the exterior unit normal to dW (see Fig. Qp. 

Notice tliat the initial-boundary value problem (|l.ip - (|1.4p is invariant under the 
self-similar scaling: 

(1.5) (t, x) — )• (oi, ax) for a / 0. 
Thus, we seek self-similar solutions with the form 



X 



for (^,77) = -. 



(1.6) (p, m, n)(i, x) = (/>, m, n)(^, rj) 

In the self-similar coordinates {(,,rj), system (jl.ip can be rewritten as 

(m - ^p)^ + {n- r]p)n + 2p = 0, 

(1.7) {p{p) - ^rn) ^ - (rim,)^ + 2m = 0, 

(Cn)g - {p{p) - 7?n) - 2n = 0. 



In the polar coordinates {r,6),r = \/^^ + the system can be further written 



as 
(1.^ 



rp — cos am — sm n ' 
dr I rm — cos 9 p{p) | + dg 
rn — s'm9p{p) 



sm am — cos n ' 
sin6p{p) 
— cos 6 p{p) 



m + '^p{p) 



n+'^p{p) 
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The location of the incident shock for large r 3> 1 is: 



(1.9) ^ = c, = JpApit:PMyo. 

V Pi- Po 

Then Problem 1 can be reformulated as a boundary value problem in an un- 
bounded domain: 

Problem 2 (Boundary Value Problem). Seek a solution of system (jl.7p . or 
equivalently (jl.Sp . with the asymptotic boundary condition when r — >■ oo; 
(1.10) 

\ (po, 0, 0) in > ,^ > 0} U {(t? - C tan 9^,)^ > 0, ?? < 0}, 

l(/3i,mi,0) in {e < 6,?? > 0}, 



{p,m,n) 



and the slip boundary condition along the wedge boundary dW : 

(1.11) {m,n) ■ u\dw = 
where v is the exterior unit normal to dW (see Fig. [2]). 

For a smooth solution U = {p,m,n) to ()1.7p . we may eliminate m and n in (jl.ip 
to obtain a second-order nonlinear equation for p: 

(1.12) ((c2 - e2)p5 - ^rjprj +^p)^+ ((c2 - r]^)p^ - ^rjp^: + r]p)^-2p = 0. 



Correspondingly, equation ()1.12p in the polar coordinates {r,9),r = Vs+^) 
takes the form 

2 2 

(1.13) {{c'-r')pr)^ + ^pr + {^pe)e = 0. 



In the self-similar coordinates, as the incident shock Sq passes through the wedge 
corner, Sq interacts with the sonic circle T sonic of state (1): r = ri, and becomes a 
transonic diffracted shock T shacks and the flow in the domain behind the shock 
and inside Tgonic becomes subsonic. In Section 2, we reduce Problem 2 for shock 
diffraction into a one-phase free boundary problem. Problem 3, for second-order 
elliptic equation in the domain with the free boundary T shock, degenerate boundary 
Tsomc; and slip boundary dWCiQ. In this paper, we focus on the existence of global 
solutions of shock diffraction and the optimal regularity of the solution across the 
sonic circle T sonic- 
There are three main difficulties to establish the global existence of solutions. The 
first is that the ellipticity degenerates at the sonic circlG Tgonic* 

The second is that 

the oblique derivative boundary condition degenerates at P2- The third difficulty is 
that the diffracted shock may coincide with the sonic circle Co := {r = c(po)} of 
state (0) in the iteration where the oblique derivative condition fails, that is, P2 may 
equal to 0. Then we can not employ directly the results in Liebermann [20]-[23] to 
show the existence of solutions for the fixed boundary value problem. One of our 
strategies here is to add an additional condition r{6) > c{po) + 6 on the diffracted 
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shock curve with 5 smah enough and modify shghtly the approximate shock curve 
to overcome the difficulty. 

The approach used in this paper for estabhshing the global existence of solutions 
is first to regularize the equation by adding the regularized differential operator eA 
to make the equation uniformly elliptic; and then to rely on the Perron method, as in 
[19] . to show the global existence of solutions for the fixed boundary value problem; 
and finally to apply the Schauder fixed point theorem to show the existence of global 
solutions for the free boundary problem. Moreover, we obtain uniform estimates for 
the global solutions with respect to 5, e > so that we can pass the limits 5 — )■ 
and e — 7- to establish the existence of solutions of the free boundary problem for 
the original system. In particular, we prove that the diffracted shock is uniformly 
transonic, that is, the strength of the shock is positive even at the point P2- 

In order to establish the optimal regularity across the sonic boundary T sonic, we 
write equation (jl.lSp in terms of the function 

:= c\pi) - c\p) 

in the (x, j/)— coordinates, which is specified in §5, defined near Pgonic sucli that Pgonic 
becomes a segment on {x = 0}, with the form 
(1.14) ^ 

(2cix — il^)tpxx + i^yy + ciV^x — '4'x ~ ~^ yj~2'^y ~ ^ in X > and near x = 0, 

plus "small" terms, since p and ■0 have the same regularity in 0,. Also -0 > in 
{x > 0} and = on {x = 0}. For our solution ip, (jl.l4p is elliptic in {x > 0}. 
It is easy to check that the solution p is Lipschitz continuous up to the sonic circle 
Tsomc = {r = ri} n dQ. 

Let W = ip — c\x. Then if) and W have the same regularity near x = 0. We 
consider the Dirichlet boundary value problem for W . We employ the approach 
in Bae-Chen-Feldman [2] to analyze the features of equation (I1.14p and prove the 
C^'"-regularity of solutions of the shock diffraction problem in the elliptic region 
up to the part rsomc\-pL of the sonic shock. As a corollary, we establish that the 
C*^'^— regularity is actually optimal across the sonic boundary rgonic from the elliptic 
region Vt to the hyperbolic region of state (1), that is, the optimal regularity at the 
degenerate elliptic boundary. 

We remark that the existence problem for a shock interaction with the right cor- 
nered wedge (90-degree) was studied by Kim [18], in which some important features 
and behavior of solutions have been exhibited. As far as we have known, for the 
shock diffraction by a convex cornered wedges whose angles are between — vr and 
in compressible fluid flow, no rigorous complete global mathematical results have 
been available, since the early work by Bargman [3], Lighthill |241 125] . Fletcher- 
Taub-Bleakney [13], and Fletcher- Weimer-Bleakney [14] . The results established in 
this paper is the first rigorous complete mathematical results through the nonlinear 
wave system for the global existence and optimal regularity of solutions of shock 
diffraction by any convex cornered wedge. 
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A closely related problem, shock reflection-diffraction by a concave cornered wedges, 
has been systematically analyzed in Chen-Feldman [6l[71[8], where the existence of 
regular shock reflection-diffraction configurations has been established up to the 
sonic wedge-angle for potential flow. Also see Canic-Keyfltz-Kim [H |5] for the un- 
steady transonic small disturbance equation and the nonlinear wave system, and 
Zheng |28] for the pressure-gradient system. 

The organization of this paper is as follows. In §2, we reformulate the shock 
diffraction problem into a free boundary problem for the nonlinear second-order 
equation (jl.ip in both the self-similar and polar coordinates, and present the state- 
ment of our main theorem for the existence and optimal regularity of the global 
solution. In §3, we first formulate the regularized approximate free boundary prob- 
lem by adding a regularized differential operator with eAp to the original equation 
(A denotes the Laplace operator in the self-similar coordinates) and the assumption 
c{p) > c(po) + ^) where p is the data given at the point P2. Then we establish the 
existence of solutions to the regularized free boundary problem for the uniformly el- 
liptic equation in the polar coordinates, and so does in the self-similar coordinates, 
as approximate solutions to the original free boundary problem. In §4, we proceed 
to the limits e — )• and 5 — )■ to establish the global existence of solutions of the 
original problem in the self-similar coordinates. In §5, we establish the optimal C*^'^- 
regularity of the solution p across the degenerate sonic boundary. In §6, we establish 
a corresponding theorem for the existence and regularity of solutions of the shock 
diffraction problem for the nonlinear wave system. 



2. Mathematical Formulation and Main Theorem 

In this section, we derive mathematical formulation of the shock diffraction prob- 
lem as a free boundary problem for a nonlinear degenerate elliptic equation of sec- 
ond order and present our main theorem of this paper. In particular, we employ the 
Rankine-Hugoniot relations to set up a boundary condition along the free boundary 
(shock) and derive other boundary conditions along the wedge boundaries in the 
polar coordinates. 



2.1. Rankine-Hugoniot Conditions and Oblique Derivative Boundary Con- 
dition on the Diffracted Shock. Consider system (jl.Sp in the polar coordinates. 
Then the Rankine-Hugoniot relations, i.e., the jump conditions, in the polar coor- 
dinates along the diffracted shock Ts^ock take the form: 

dr f sin 6* cos 9 \ 

(2.1) —f{P\ + cos t'[mj -|- sm.d[n\ — '^y "I M j ; 

_r -, , , dr sin 6 

(2.2) cose/ p - r[m\ = - — p , 

dO r 

■ 1 r 1 dr cos 9 , , 

(2.3) sm9[p]-r[n] = — \p]. 
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Solving for [m] in (j2.2p and for [n] in (|2.3|) . we have 
, , , . , /dr sin 9 cos 6 \ , , . , / d 

H = (5^— + — )H. ["1 = (-;j 

Combining (12. ip with (I2.4p . we find that, for our case, 



with c{p, po) = xi ^^^pJ'p^^^ , where we have chosen the plus branch so that % > 0. 



Moreover, (I2J)-([MD imply that 

(2.6) \p][p] = H' + [n?. 
Differentiating (j2.6p along Tshock) i-e-, r'{9)dr + 80, yields 

(2.7) {c^[p] + [p]) (p,r'(0) + pe) = 2H(m,r'(0) + mg) + 2[n](ny(0) + ng). 
From the second and third equations of (jl.ip and uj = — w-^j , we obtain 

ut = 0, 

which means that, in any region formed by the flow trajectories across a curve along 
which the irrotational condition oj = holds, the classical solutions must be also 
irrotational. That is, = holds in the region up to the shocks in the self-similar 
coordinates. Then, in such a region, using (|1.8p . we have 



(2 



2 2 

mr = ^ cos 9 pr- yi sin 9 Pe, 

2 2 

nr = ^ sin 9 pr + ^cos9 Pe, 



me = (c^ — r^) sin 9 p^ + ^ cos 9 pe, 

ne = -(c^ - r"^) cos 9 pr + ^sm9pe. 

Hence, from (j2.7p . replacing {mr,me) and {nr,ne) by the above equations and using 
[m] and [n] in (|2.4p . we obtain 

2 

(2.9) (3iDip := + /Sa/Og = 0, 

1=1 

where /3 = {Pi, 132) is a function of {po, p,r{9),r' (9)) with 

/3i = r'(0)(c2(r2 - c2) - 3c2(c2 - r^)), 

/32 = 3c2(r2-c2)-c2(c2-r2). 

Thus, the obliqueness becomes 

/ /3 . (1, -r'(0)) = -2r\c^ - c'^)r'{9) =: p, 

where (1, —r'{9)) is the outward normal to Q on Tshock- Note that p becomes zero 
when r'{9) = 0, that is, r = c{p,pq). When the obliqueness fails, we have 

/3i = 0, f32 = -c\c^ - r^) < 0, 



8 GUI-QIANG G. CHEN XUEMEI DENG WEI XIANG 

since c^{p) > c^{p,po) = r"^ if p > po. 

We define Q to be the governing second-order quasilinear operator in the subsonic 
domain 0: 

2 2 

(2.11) Qp := ((c2 - r^)p,)^ + ^pr + i^pe), = 0, 
and M to be the oblique derivative boundary operator: 

(2.12) Mp := Pipr + hpe = on T.hock := {(r(^), ^) : 9.^ < 9 < 
The second condition on Fshock is the shock evolution equation: 

(2.13) % = ,V^EF^^ ^(,, 0^ ^))^ ,(0^) = 
d9 c{p,po) 

where (ri,0i) are the polar coordinates of Pi = {Ci^Vi)- 

At the point P2, r'{9uj) = 0, M does not satisfy the oblique derivative boundary 
condition at this point. We may alternatively express this as a one-point Dirichlet 
condition by solving r{9u}) = c{p{r{9w), 9^), po)- In order to deal with this equation, 
we introduce the notation: 

(2.14) a = (cb)~^(r) when Cf, := c{a,b) = r for fixed b. 
Thus, we have 

(2.15) p{P,) = p={cp,)-Hr{9^)). 

2.2. Boundary Condition on the Wedge. The boundary condition on the wedge 
is the slip boundary condition, i.e., 

(m, n) ■ u = 0. 

In order to derive the boundary condition for p itself, we notice that the normal 
direction along the wedge is (— ™d the tangential direction is (^,?7), which 
implies that — rjm = 0. Differentiating it along the wedge yields 

= ^^ng — ^rjm^ + + ^rjUr^ — rf'nir^ — rjm = ^'^n^ — ^rjm^ + ^fl^r] — fffriri. 

Combining this with the second and third equations in (jl.ip . we conclude that p 
satisfies 

(2.16) pu = ^ onVQ:= dVt.r\ {{9 = -k}\J {9 = 9 u,}). 

2.3. Boundary Condition on ^ sonic State (!)• The Dirichelt boundary con- 
dition on T sonic- 

(2.17) p = pi on Tsonic := n dB,M- 

On the Dirichlet boundary Fsonio the equation Qp = becomes degenerate elliptic 
from the inside of 
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2.4. Reformulation of the Shock Diffraction Problem. With the derivation of 
the free boundary condition on Tshock and the fixed boundary conditions on Tgonic 
and the wedge To, Problem 2 is reduced to the following free boundary problem 
in the domain Q for the second order equation (j2.1ip . with {m,n) correspondingly 
determined by (12. 8j) . 

Problem 3 (Free Boundary Problem). 5*66^ a solution of the second-order 
nonlinear equation (|2.1ip for the density function p in the domain ft, satisfying the 
free boundary conditions (I2.12p -( r2.15p on Tshock, the Neumann boundary condition 
(j2.16p on the wedge Tq, and the Dirichlet boundary condition (|2.17p on the degenerate 
boundary Tsonic {the sonic circle of state (1)) (see Fig. [2]). 

2.5. Main Theorem. For the free boundary problem, Problem 3, we have the 
following results, which form the main theorem of this paper. 

Theorem 2.1. (Main Theorem) Let the wedge angle 9w be between — vr and 0. 
Then there exists a global solution p{r, 9) in the domain Q with the free boundary 
r = r{6),ee [eni,Oi], o/ Problem 3." 

p e C72+"(17) n C7"(n), r G C2+"([0^, 9,)) n C''H[9^,9i]). 

Moreover, the solution {p{r,9),r{9)) satisfies the following properties: 

(i) p > pq on the shock Tghock, that is, the shock Tghock is separated from the 
sonic circle Cq of state (0); 

(ii) The shock T shock is strictly convex, except the point Pi, in the self-similar 
coordinates (^, r]); 

(iii) The solution is C^'°' up to Tsonic cmd Lipschitz continuous across Tsonic^ 

(iv) The Lipschitz regularity of solutions across Tsonic o,nd at Pi from the inside 
is optimal. 

In particular, Theorem 12.11 implies the following facts: 

(i) The diffracted shock T shock definitely is not degenerate at the point P2. This 
has been an open question even when the wedge angle is ^ as in [18], though 
it is physically plausible. 

(ii) The curvature of the diffracted shock Tshock away from the point P2 is strictly 
convex, though the strict convexity of the curvature fails at P2 . 

(iii) The optimal regularity of solutions across Tsonic and at Pi from the inside 
is C^'^, i.e., Lipschitz continuity. 

We establish Theorem 12.11 in two main steps. First, we solve the regularized 
approximate free boundary problem for Q involving two small parameters e and 6, 
introduced in §3. Then we analyze the limits e — t- and 5 — t- 0, and prove that the 
limits yield a solution of Problem 3, i.e., ()2.1ip - ()2.17p . in §4. The optimal regularity 
is established in §5. 
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3. Regularized Approximate Problem 

In this section we first formulate the regularized approximate free boundary prob- 
lem and establish the existence of solutions to this problem as approximate solutions 
to the original problem. 

To solve the free boundary problem, we formulate the fixed point argument in 
terms of the position of the free boundary. There are two main difficulties in es- 
tablishing the existence of solutions: The first is that the ellipticity degenerates at 
the sonic circle rgonic 

; and the second is that the free boundary Fshock may coin- 
cide with the sonic circle Cq of state (0) during iteration, which makes the iteration 
impossible. We overcome these difficulties as described below. 

3.1. Approximate Problem and Existence Theorem for Approximate So- 
lutions. For fixed e, introduce a regularized operator 

Q' := Q + eA, 

where A represents the Laplace operator in the self-similar coordinates. For a given 
curve r{9), we first solve the fixed boundary value problem (j2.1ip - (|2.12p . (|2.16p - 
(j2.17p . and (j2.15p with replacing Q; then we obtain a new shock position r{6) by 
integrating (|2.13p : 

(3.1) r{e)=n+ [ g{r{s),s,pis,ris)))ds for 9 e [6^,61), 

Jei 

where g is defined in ()2.13p . Note that, on the right side of (|3.ip . we evaluate all 
the quantities along the old shock position r{6). 

With this, it seems that the free boundary could be obtained by solving a fixed 
boundary problem and then by integrating the shock evolution equation. However, 
we face the second difficulty as indicated above, that is, f{9) may meet the sonic 
circle Cq of state (0). Introduce another small, positive parameter 6 which is fixed 
and define the iteration set of r, /C^'^, which is a closed, convex subset of a Holder 
space C^^°''^ {[6w,0i]), where ai depends on e and 6 to be specified later. The 
functions in /C^'^ satisfy 

(Ki) r{ei) = n- 
{K2) r'(^^) = 0; 
(K3) c{pQ) + 5<r{e^)- 

(Ki) 0<r'{e)<£^ iov6^<6<9i. 

When the difficulty occurs, we modify r{6) slightly somewhere as r{6) = c(po) + S + 
A{6 — 0«,)^ + B{6 — 0«,)", where A, B, and n will be uniquely determined. Then we 
define a mapping on JC^'^: 
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We now restate the regularized approximate problem as follows: For fixed £,6 > 0, 
the equation for p in the subsonic region is 

2 2 

(3.2) Q^p = ((c2 - + e)pr)^ + + i^^Pe)e = 0; 

the shock evolution equation remains the same when r > c{pq) + 26: 

(3.3) 



-9{r,9,p), 



) de 

\r{Oi) = ri; 
and 

(3.4) r{e) = c{po) + 6 + A{d- O^f + B{e - O^T 

for some constants A,B, and n on the boundary when (13. 3p does not hold; the 
remaining boundary conditions as before are 

(3.5) Mp = /3 • Vp = on Tshock = {(r, 9) : < < ^i}, 

(3.6) P = Pi on Tsonic; = on Tq, 
where v is the outward normal to at Tq; and 

(3.7) p{P^)=-p={cp,r\r[e^)). 

_ Let 1/ = {Pi,P2,0,i^3} denote the corners of 0, and V = V\{P2}. Set Q' = 
0\(y U Tshock)- For P G y, we define the corner region 

VLp{(t) := {x eVL : dist(x,P) < a}, ily(cr) := Upeyl7p(cr). 

We define a region that is close to Fshock) but does not contain the corner Pi by 
taking a covering of Fghock with a ball of radius 5 centered at the points on Fshock 
which are bounded away from Pi. Define 

T'ia) := {P G Fshock : dist(P,Pi) > a} 

and 

T{a) = {xenn{Up^r'i.) BAP))}, 
where B^{P) is a ball of radius a centered at P. We then define 

(3.8) = {u : \\uf^ := sup {c7''+^\\u\\afl\(r{a)uny,{a))) < oo}. 

cr>0 

We focus now on the proof of the following existence theorem in this section. 

Theorem 3.1. For any e G (0, eo) f^n-d S G (0,5o) for some £o,5q > 0, there exists 
a solution (p^'^, r*^'"^) G C'^^"^^{^}'^'^) x C^+°([0^, ^i]) to the regularized free boundary 
problem ()3.2p - (j3.7p such that 



(3.9) Po<p'''^<p'''^<Pi, c^ip''^)>r'^ mn''^ 

for some a, 7 G (0, 1), which depend on e, 6, and the data {po, pi,Ow)- Furthermore, 
the solution satisfies (j3.3p at the points ofT^^f^^^j^ where r"^'^ > c(po) + 26. The curve 
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r^' (9), defining the position of the free boundary T^^^^^, is in JC''^. Here Q''^ IS 
bounded by r^jf^cfc' ^ sonic, and Tq. 

We establish Theorem 13.11 in the following steps whose details are given in the 
following four subsections. 



Step 1. Since the governing equation (j3.2|) is nonlinear and the ellipticity is not 
known a priori, we impose a cut-off function in the equation p = 0. 
We introduce a smooth increasing function Q G C°° such that 



(3.10) as) :-- 



s if s > 0, 
— ^e if s < — e, 

and ^'(s)! < 1. We then consider the following modified governing equation: 

Qe,+p = ^^^^c^_r^) + e)pr)^ + i^pe)9+{m^'-^') + ^) + ^)Pr 

(3.11) 

= E'i=iH''Ur,e,p)D^p) +b'{r,p)DrP = in fi. 

Step 2. We make some estimates for a solution to the linear problem with fixed 
boundary Fshock defined by r{0) G /C^'^ and establish the Schauder estimates on 

Tshock' 

For the linear problem, we establish a priori Schauder and Holder bounds on 
Tshock, especially near the point where the boundary condition loses the obliqueness. 
We use the Holder estimates of the neighborhood of the corners and the C2+"- 
estimates locally in the rest of the domain. 



Step 3. We employ a technique in Lieberman |19j to solve the problem with the 
oblique derivative boundary condition Mp = 0. Using the Holder gradient bounds 
to the linear problem, we establish the existence results for the linear fixed boundary 
problem in the polar coordinates, via the Perron method developed in |19j . 



Step 4- We apply the Schauder fixed point theorem to conclude the existence of 
a solution to the free boundary problem with the oblique derivative boundary con- 
dition. Finally we remove the cut-off functions by the a priori estimate to conclude 
the results. 

3.2. Proof of Theorme 3.1: Regularized Linear Fixed Boundary Value 
Problem. Replace p in the coefficients af^ and b^ in ()3.1ip and /Jj in (13. 5p by a 
function w in a set W that is defined with respect to a given boundary component 
■'^shock depends on the given values po, pi, and p^'^ = {cpg)~^ {r^'^ (9^)): 

Definition 3.2. The elements w oiW C C? ^ satisfy 

(-71) 

(Wl) Po < p"^'^ < w < pi, w = pi on Fsonic, = on Fq, and w{P2) = p^'^; 
(W2) \\w\\ao < Ko, ||u;||2+aoAoc ^ ^'0, and \\w\\i+^^rido) < Kq] 
(W2) \\w\\i~^''' < Ki. 
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The weighted Holder space is defined by (j3.8p . The values of 71,00, /i G (0,1), 
and do, as well as the values of Kq and Ki, will be specified later. Obviously, W is 
closed, bounded, and convex. 

The quasilinear equation (j3.1ip and the oblique derivative boundary condition 
(|3.5p are now replaced by the linear equation and linear oblique derivative boundary 
condition on T^J^^^^ := {{r{9),9) : 6^ < 6 < 9^}: 

L^^+u := Y:Ui Di{al-{P, w)Diu) + ¥{P, w)Diu = in 

Mu:= /3i{P,w)DrU + /32{P,w)Deu = on Tl^^^^, 

with given r{9) G IC'^^ C ([0^, n C^{{9^,9i)) and w £ W, where the 

repeated indices are summed as usual. Because of the cut-off function L^'~^ is 
uniformly elliptic in ft with the ellipticity ratio depending on the data and e. 

In this section, we demonstrate the key point that, for a given function w G W, 
the solution u to the linear equation (|3.12p with the remaining boundary conditions: 

(3.13) u = pi on Tshock, Uj, = on To, u{P2) = p"'^, 

satisfies the Holder and Schauder estimates in Q,' and the uniform bound in C^~^'^(r{dQ)) 
near T^^^^^j^ for any fi < min{7i,ai}. This bound gives rise to enough compactness 
to establish the existence of a solution to the quasilinear problem by applying the 
Schauder fixed point theorem. 

First, we state the Schauder estimates up to the fixed boundary Tsonic with the 
Dirichlet boundary condition, to To with the Neumann boundary condition, and the 
Holder estimates at the corners V' . 

Lemma 3.3. Assume that Tghock is parameterized as {(r{9),9)} with r{9) G /C"'^ 
for some ai and that w G W for given Kq, Ki, cko, and 7. Then there exist 
IViOin G (0,1) such that any solution u G Cf^t""(^') ^ C^^ {^v'{d(i)) to the linear 
problem (I3.12p - (l3.13p satisfies 

(3-14) \\u\\'r,ny,{do) < Ci\\u\\o for any 7 < 7y, 

and 

(3.15) ll«ll2+a,Q; <C'2||ii||o for any a < an. 

The exponent depends on the data pQ, pi, and 9^; and both an and 7v depend 
on £ but are independent of ai and 71 . The constant C2 is independent of Ki but 
depends on Kq. 

Proof. The corner estimates at Pi and P3 directly follow from the results in Theorem 
1, Lieberman [23]. Near the origin, the wedge angle is larger than vr; thanks to the 
symmetry of the governing equation in the 0-axis with form (jl.l3p . we derive the 
corner estimate as follows. 

First we fiat out the boundary by introducing the transformation: 

{r',9') = (r,^— (0-0^)), (e',7/) = (r'cos0',/sin0O. 
vr — U,n 
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Then the governing equation in the (r', 0')-coordinates takes the form 



-'w 

and 

+ (((C(cV) - + - 

in the (^', ?7')-coordinates, where we drop ' for simphcity without confusion. The 
eigenvalues of ()3.16p are 

Ai = C(c' - r') + e, X2 = {^^)\c' + e). 

If we can show the transformation from the r/)— coordinates to (^', 77')— coordinates 
is invertible and the C°-norms are equivalent, then we can obtain the estimate of 
the solution at O by applying the standard Schauder-estimate since the coefficients 
of Q^'^ p are in L°° . In fact, 

^ = cos COS (-^ (0 + 0o)) + sin sin (-^(^ + ^o)) , 
% = si^^^os {^^(^ + ^0)) - cos sin (^(^ + 0o)), 
% = ^o,9 sin (^(0 + ^0)) - sin cos (^^(^ + ^o)) , 
% = ^ sin + ^0)) + ^ cos cos (^(^ + 0o)) . 



Thus, 



which implies that the transformation between the two coordinates is invertible. 
As for the proof of the equivalence of the two norms, we have two cases: 

Case 1. If 6* > f as in Fig El then 



e' = ke>- with k ^" 



2 TT 

Since 1 < /c < 2, 

\x — y\> max{r(x), r(?/)}, \x' — y'\ > max{r(j;'), r(y')}. 

Then the equivalence of the two C^'^— norms can be easily shown by setting r{x) 
r(x') and r{y) = r{y'). 



SHOCK DIFFRACTION BY CONVEX CORNERED WEDGES 



15 




Figure 3. Scaling of the angles 



Case 2. li 9 < ^, then the distance between two points in the ry)— coordinates 
and (i^', ry')— coordinates is equivalent. By the cosine law, we define 

m{k) := \x' — y'p = r(x)^ + r(y)^ — 2r{x)r{y) cos{k9), 



and then 



Thus, 



^^^^ = 2kr{x)r{y) sm{k9) > 0. 

OK 



\x — y\ = \J m(l) < \J m{k) = \x' — y'\ < \J m{2) < 2\J m(\) = 2\x — y\. 

Therefore, we can obtain the Holder estimate of the solution at O. Here 71/ de- 
pends on the angle at the corner, a fixed value that depends on the data {po, pi,9w), 
and the ellipticity ratio e, but independent of 71, ai, Kq, and Ki. 

Finally, we can use the standard interior and boundary Schauder estimates to 
obtain the local estimate ()3.15p . The constant C2 depends on e, the C"— norm of 
the coefficients Qjj, and the domain. □ 

Because the interior Schauder estimates can be further applied, a solution in 
Cl+'^iQ') is actually in Cl^{n). 

We next establish the Holder gradient estimates on Fghock- It is at this point that 
we need to derive the basic estimates at the point P2 where the boundary operator 
M is not oblique. In order to avoid handling the Neumann boundary condition on 
the wedge boundary 9 = 9w separately at each step of this proof, we reflect $7 across 
the wedge boundary 9 = 9^, without further comment, i.e., O includes Sq, and 
let 

Fshock stand for the full C^^*^^— boundary in Lemma |3i4j below. In addition, we 
extend u{29w — 9) = u{9) for 9 G {9w,9i) in a small neighborhood of 9w We still 
denote n by u for simplicity without confusion. 

Lemma 3.4. Assume that Tghock is given by {(r(9),9)} with r{9) G /C^'^ for some 
«! and that w G W for given Kq, Ki, oq, and 71. Then there exists a positive 
constant do such that, for any d < do, the solution u G C\^^{U, U T shock) U (^^^^(O) to 
the linear problem (I3.12p -( r3.13p satisfies 

(3-17) lklli+M,r{d)\Bd(Pi) <C{e,6,ai,-^i,Ki,do)\\u\\o 
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for any fi < min{7i, ai}. 

The proof of this lemma is adopted from which is long and tedious, so we 
postpone it in Appendix for self-contained. 

Now we focus on the existence of solutions in Theorem 3.1 for problem (13.12p - 

First we introduce two definitions with some modification in comparison with |19j . 

We say that problem (j3.12p - (j3.13p is locally solvable if, for each y £ il, there exists 
a neighborhood 0{y) such that, for any h G C(N) with N := 0{y) n {f7\({P2} U 
rsonic)}) there exists a solution v € C'^{N) n C{N) of the problem 

L^^+v = in n 0, 
Mv = Q on iV n 90, 

V = h on d'N, 

when P2 ^ N{y); or 

L^'+u = in A n 17, 
Mv = on A n dn, 

V = h on d'N, 
v\p, = P''', 

when P2 £ A"(y). Here d'N = dN D O. For brevity, we denote this function v by 
{h)y to emphasize its dependence on h and y. 

A subsolution (supersolution) of (|3.12p - (j3.13p is a function w G C(J7) with 
'w{r{9ii;),9uj) < /J^'*^ {w{r{6u)),9u)) > p'^'^) such that, for any y G $7, if /i > u; 
{h < w) on d'N, then (/i)^^ > w {{h)y < w) in A. The set of all subsolutions 
(supersolutions) is denoted by (5*^). 

We now establish the existence of solutions to problem (j3.12p - (j3.13p . 

Lemma 3.5. Assume that Tghock is given by {{r{9),d)} with r(9) G /C^''^ for some 
ai and that w gW for given Kq, Ki, uq, and 71. Then there exist 7\/,an G (0, 1) 
and do > 0, which are independent of and ai, such that there exists a solution 

to the linear problem (j3.12p - (j3.13p for any a < an, p < min{7i,Q!i}, 7 < 7^, and 
d < do, which satisfies (j3.14p - (|3.15p and (|3.17p . 

Proof. For fixed e,S > 0, we denote u"^'*^ = u in the proof without confusion. We use 
the Perron method to show the existence of a solution to problem (j3.12p - (j3.13p . 

It suffices to show the local existence at P2. In fact, let B2 be a sufficiently small 
neighborhood of P2 with smooth boundary such that O ^ B2, Pi < 0, and (32 < 0. 
Then we study the local existence in the (,^, r/)— coordinates in B2- We introduce 
the coordinate transform in a neighborhood of P2: 

(3.18) i = i{r,e), f, = f,{r,e) 
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such that 

Cirw, Ow) = 0, r/(r^, 9w) = 0, 

^ = -1 ^ = > n 

dr ' 9^ ^2 ~ ■ 

Let Tshock := {{Lv) : f, = /(!)} = {(r,0) : r = r(0)} in B2. Then r?(r(e),0) = 
f{i{r{e),9)) and hence 

and the function /(^) is increasing in | on Tshock H -B2. Thus, from || = — -g^ > 
and 1^ = 0, we have 

fii) > 0. 

We reflect the region B2 across ^ = to obtain a new region, which is still denoted 
by B2. Furthermore, we replace by which is the a— distance from the point 
P2 upward, see Figure HI On the bottom straight boundary of Og-, we impose 

u = p^'^ on bottom of 0,^- 




1 




Figure 4. Domain with tip P2 removed 
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Now we study the following boundary value problem: 

' L'^'^u = YA,j=i aijDjU + Ylf^i biDiU = inil^, 



(3.19) 



Mu = d^U = on dVLa H Fshock! 

u = h ondB2 n O, 

u = p^'^ on So-) 



where 



«11 - ^) «12 - «21 - -3^a|2, 022 - Oil + (y-)^a22, 

P2 P2 P2 



' /3i 5e /3i /32 di PI df, ■ 

Here afj, 6"^, and z = 1, 2, are the coefficients of (|3.12|) - (|3.13|) in the ?})— coordinates, 
and h is a continuous function satisfying p^'*^ < h < pi. Following Lieberman |19j . 
there exists a solution 

Mo G C(nn S2) U C72'"(0^ n S2) 

for B2 small enough. The maximum principle holds for which converges locally 
in C^{n n B2) to a solution in C^+'^in D B2) as cr ^ 0+. 

We now use a barrier function to obtain the continuity of u at P2- We consider 
the auxiliary function 

(3.20) v = p^^^ + c{l-e-^^), 

where c > and / > are specified later. For the oblique derivative boundary 
condition along n Fghock, we have the following two cases: 

Case 1: ■ v > and M{v - p^'^) > when | > 0; 
Case 2: ■ u < and M{v - p"'^) < when | < 0, 

where v denotes the outward normal to Oo- at ^o- H Fghock- 

Moreover, it is easy to see that v > p^'^ on So-. Choose C large enough such that 
V > sup \ h\ on dB2 Ci For the equation, we have 

dIjDjV + bIDiV = -ce-^^{fh\^ -ib\), 

which is less than a negative constant by choosing / > ||6^||o/A, where A < af]^(^, ?]). 
Thus, 

p"'^ <u<v. 

Then u is continuous at the point P2- The continuity of u at the other points follows 
from Lieberman's argument in [211 123j . 
By Lemma [Ql we have u £ C^''^{B2 nH). 
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In order to establish the global existence of solutions, it is required to show 

sup(tt;~ — w'^) = 0, 
n 

where are the supersolution and subsolution of problem (j3.12p - (j3.13p . respec- 
tively. 

In fact, we set m := sup(it;~ — w~^). We assume that m > in 0. Since w~{P2) — 
n 

w~^{P2) < 0, there exists a neighborhood B2{P2) of P2 such that (y) — (y) < m 
for y G B2{P2)- Now we define 

3^ := {y G O : w~ (y) — w^{y) = m}. 

Let yo £ y such that 

dist(?/o,-P2) = mindist(y,_P2)- 

Let w"^ be the lifts of w"^ in M{yo). We see that 

iv~ — < m on d'N. 

The strong maximum principle implies that either w~ — w~^ < min M or —w'^ = 
m. Since 'w~{yo) — w'^{yo) > w~ (yo) — w^{yo) = m, it follows that 

w~ — = m in A^, 

and hence 

iB~ — w'^ = m on 9'iV, 
which contains the point of y closer to P2 than yo- This is a contradiction with the 
definition of y^. 

We refer to Lieberman [21] to handle the mixed case and the points Pi and P3, and 
Lieberman |22j to handle the point O where the two Neumann boundary conditions 
are satisfied. As for the interior and the Dirichlet boundary condition on the sonic 
arc Vsonic^ they are classical since the equation is uniformly elliptic for fixed e > 
(see Gilbarg-Trudiger [16]). 

With all of these, we then employ the Perron method to establish the existence 
of a global solution. □ 

3.3. Proof of Theorem 3.1: Regularized Nonlinear Fixed Boundary Prob- 
lem. We now establish the existence of solutions to the nonlinear problem ()3.2p 
with a fixed boundary. 

Lemma 3.6. For e G (0,eo) and 5 G (0,(5o), given r{9) G /C^'^ C C^+"S there 
exists a solution p^'^ G C'^^"^s^{Q^'^) to problem (j3.2p and ()3.6p - (|3.7p with the oblique 

derivative condition Mp^'^ = for some a(e, 5), 71 (e, 5) G (0,1) such that 

(3.21) Po < p^'^ < p''^ < pi. 

Moreover, for some do > 0, the solution p^'^ satisfies 

(3-22) \\p''^\U,r{do)uBag(Pi) <-^2, 

where 7 and K2 depend on 5, e, ^y, and Ki, but are independent of ai. 
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Proof. For the notational simplicity, we write p = p^'^ throughout the proof. 
For any function E W, we define a mapping 

by Tw = p, where p is the solution to the linear regularized fixed boundary problem 
()3.12p - ()3.13p solved in Lemma [3.51 Because of the cut-off function L^'+ is strictly 
elliptic, with the ellipticity ratio depending on e. By Lemma |3.5| T maps W into 
a bounded set in C'^'^^^y where 7y is the value given by Lemma [331 Since 7^ is 

independent of 71, we may take 71 = 7v/2 so that r(/C) is precompact in C'^_^_^y 

To show that T maps W into itself, we need to show that Tw satisfies (Wl), 
(W2), and (W3) in Definition O 

Now, (Wl) is immediate by the boundary conditions and the maximum principle. 

By applying the standard interior and boundary Holder estimates (cf. Theorems 
8.22 and 8.27 in [K]), we obtain 

(3.23) \\p\L*,n', < Co, 

where a* G (0, 1) and Cq depend only on e (the ellipticity ratio), d' = dist(r2'^, 50') 
with Q'l C ri', and the data po, pi, and 6w Notice that the constant Cq is non- 
decreasing and the constant a* is nonincreasing with respect to d' (see the remark 
following Theorem 8.24 in [16|). Since fi' C is bounded, we can find an upper 
bound for Cq and a lower bound for a* depending only on the size of and the 
ellipticity ratio. Thus, if we define W with Kq = Cq and ao = a*, then p = Tw 
satisfies (W3). Note that Kq and a* are independent of ai and 'jy- 
To verify (W2), it suffices to find a constant K such that 

(3-24) sup (52-71 ||p||2^f^\|p(5)un„.(5)}) < K, 

under the assumption that \\w\\2 "'^'^ < K. Note that Lemma 13.41 implies the exis- 
tence of a constant do > such that, for d < do, any solution u £ C^{il:U Fghock) U 
C^(0) to the linear problem (j3.12p - (j3.13p satisfies (|3.17p . where the constant C de- 
pends on K but is uniform in d < do. Based on this estimate, we obtain a local 
bound for the weighted norm of p on T{dQ) of the form 

(3.25) d^-^'\\p\\2<d^-^''+^'C, 

which holds for all d < do, where C depends on K, ai, and 71. To show (|3.24p . we 
make the L°°-estimate by considering separately the domains in r2\{r(5) U Qv{^)} 
for which 6 > d, with d < do to be specified later, and the domains for which 5 < d. 

In the domains of the first kind, 0\{r((5) U ily((5)} with 6 > d, the solution is 
smooth, and its C^-norm bound is independent of K. More precisely, we can use 
the uniform Holder estimate (|3.23p and the bootstrap iteratively to obtain 

(3.26) \\p\\2+aa,n' < C{Ko). 

Notice that, since the Holder estimate ()3.23p is independent of the distance between 
^'i and the boundary Fshock, so is the Schauder estimate ()3.26p . The interpolation 
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inequality (cf. Lemma 6.32, [16]) gives 

(3-27) \\ph,n' < c\\p\\o + cr\\p\\2+an,n' < cpi + aC{Ko) 

for any a > and c = c{a). We fix o" = 1 and find 
(3-28) sup (52-71 ||p||2,n\{rwuQH5)}) ^ 

5>d 

where K' depends on the size of the domain Q, C{Ko), and pi, but is independent 
of the distance to Fshock- 

Next we estimate 6 '^^ IIPII2 f2\{r((5)uf2v'('5)} ^^^^ 6 ^ d. We divide the subdomain 
r2\{r((5)uOy((5)} into two parts: The part for which 6 > d and its complement. The 
supremum over the subdomain for which 5 > d has been calculated above. We use 
the estimates for the behavior of the solution near Fshock to obtain the supremum 
over the complement, that is, estimate (j3.25p and the corner estimate (j3.14p . In 
(|3.14p . the constants Ci and 7y are independent of K, Kq, and ai, where \\p\\o is 
bounded by pi. By the interpolation inequality, we have 

(3-29) WpW^yuUvidv) ^ CvWpW^y^nvidv) ^ CyCipi, 

since 71 = 7\//2, where Cy = Cv{'yi,^v,^vidv)) for some dy > 0. From here, we 
obtain 

d'^~'"\\p\\2<Kv, yd<dv, 

where Ky is independent of K. Let K = max{Kv,K'}. Using bound (|3.28p . K 
is independent of qi and d. Since Ky and K' are independent of d, we change d 
without affecting K. Therefore, we can choose d < minjdo, d\/}/2 in ()3.25p small 
enough that 

Jl-7i+M(7 < K. 

Therefore, (j3.24p is satisfied, and we have chosen the parameters K, Kq, and oq 
defining W so that T maps W into itself. 

Now, by the Schauder fixed point theorem, there exists a fixed point p such that 
Tp = p G C'^_ryj^y Thus, /3 is a solution of the boundary value problem (13. 2p and 
(|3.5p - (j3.7p . By a bootstrap argument, we obtain 

p e fo'^ « < «n, 

where ai is given in Lemma 13.31 Note that we have chosen 71 = 7y /2 for the 
exponent 7y € (0, 1) depending on e and the corner angle an. 

The maximum principle gives the first estimate in (j3.2ip . and the second follows 
from Lemma [ 



Finally, since T(W) C W is a bounded set in C'^_^_^y by the interpolation in- 
equality, any fixed point p satisfies (j3.22p for any 7 < 71. Note that K2 and 71 are 
independent of ai. □ 
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3.4. Proof of Theorem 3.1: Three Important Properties for the Nonhn- 
ear Problem. In this subsection, we show three properties for the fixed boundary 
nonhnear problem (13. 2p and (I3.5p -( f37f|l . We first prove the following inequality: 

(3.30) c\p''^)-r^>0 mn"'\ 

which guarantees the ellipticity of our nonlinear system, so that we can remove the 
cut-off function. 

Lemma 3.7. There exist positive constants Eq and Sq such that, for any fixed < 
e<£o andO<6< 5o, the solution p^'^ G CmnC'^{n)r\C^{n\T sonic) to (IO)- (f3Tll 
satisfies 

(3.31) c^{p''^)>r^ inQ''\ 

Proof. For the notational simplicity, throughout the proof, we write p = p'''^. 

We prove the result by contradiction arguments. On contrary, we assume that 
there exists a nonempty set Z) = {(^, 77) G : c^{p) — < 0}, and choose Xynm £ D 
to be the minimum point. Then it is easy to check that P2 ^ D. Since O ^ D, 

Dcns:={X£ n\V : r^ > c\p,po)}, 

where V is the set of the corner points of ^l. Thus, there are three possible locations 
of X^[ji. 

Case 1: Xmin G D locates inside Q,s- For notational simplicity, we denote c^(yo) = 
p'*'"^ = u. Then, multiplying (7 — l)p'^~'^ both sides of the equation Q^'^p = 0, we 
have 
(3.32) 

Lu = (7 - l)p^-2ge,+^ 

= Ell aUDnU - Anp) + C'(c2 - r^)ic' - r\ur + ^uj + ¥ur 

= 0. 

We note that | < ofi < £ due to the cut-off function ( in D. We evaluate Lr^ in D: 

Lr^ > -2e\l-^^-^r^\+C{c^-r^)ic^-r\ur + 2c^ 
^^'^^^ > 2r2-^|n^-l-fc^r2| >0 

Pq I 

with small e < Eq, where 

7-1 2 
Po ^0 



17-1 _ 2|7-2| 21 
|/^0 7-1 ^01 



and the fact that (c^ — r^)r.(Xniin) = 0. Then, in D n $7, we obtain 
(3.34) 

> Lu - Lr^ 

= E-=i {<^Du{u - r2) - ^_I^af,A(n + r^)D,{u - r^) 

+C'{c^ - r2)(c2 - r%{u - r'^)r + X{u + r'^)g{u - r'^)g + b{u - r% 



SHOCK DIFFRACTION BY CONVEX CORNERED WEDGES 23 
Since Xmin is an interior minimum point, 

2 

A {U - r^){X^in) =0, Yl - (^min) > 0, 

1=1 

which contradicts the inequahty Lu — Lr"^ < in D n 17. 

Case 2: Xmin locates on Fshock H D. Multiplying (7 — l)p^~'^ over the equation 
Mp = 0, we have 

2 

= (7 - l)p^-'^Mp = Mu = Y (3iDiU. 

i=l 

On one hand, we have 

(3.35) Mr^ = 2r/3i = 2rr'{c^{r^ - c^) - Sd^ic^ - r^)) > on r^hock n D, 

where we have used the fact that > > (P in Qg. 

On the other hand, at Xyaim the outward normal derivative of u — becomes 
non-positive (i.e., V(u — r^)(l, — r') < 0) and the tangential derivative becomes zero 
(i.e., V{u - r2)(r', 1) = 0). We obtain 

{l + {r'f){u-r\<0 atX^in, 

which implies (n — r^)^ < 0. Thus we have 

> M{u - r^) = (/3i - r'f32){u - r^)r = p{u - r^)^ > 0, 

which is a contradiction. 

Case 3: Xmin locates on TqCi D. Then 

which is a contradiction due to the Hopf maximum principle, i.e.. 

Therefore, there is no minimum point, which implies the set D = 0. This com- 
pletes the proof. □ 

Based on Lemma [3. 71 we can show that the solution of the free boundary problem 
(I32])~(I3Z1) satisfies 

r - c{p,po) > on Fshock- 

Lemma 3.8. Let < £ < £0 and < 5 < 6q, and let 

p''^ G c{n) n c\n) n c\n\Tsomc) 

be a solution of the boundary value problem (13. 2p and (I3.5p - p.7p . Then 

c{p''\po)-r<0 onTli^,. 
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The basic idea of the proof of this lemma is from ^E\. For self-containedness, we 
include the proof in Appendix. 

Based on Lemma [3^ it is easy to show that r{9) > c{p'^'^ {r{9),6), pq) on T^^^^j^ 
for 6 G [Oyj^Oi). Thus, we have shown that the integration in (jS.ip can always be 
integrated on Tshock- 

In the next lemma, we study the monotonicity of p along the shock boundary 
Tshock) which will be used to describe the behavior of p^'^ and r^'*^ near the shock 
-'^shock ■^hen e, 5 tend to zero, and the convexity of the shock in the r/)-coordinates. 

Lemma 3.9. Suppose that p^'"^ G {Q.[JT shock^'^ wedge^^Q)^C"- {^) is a solution of 
the boundary value problem (13. 2p and ()3.5p - ()3.7p . Then p""'^ is monotone on T shock- 

The idea of the proof of this lemma is from [5], and the main difference here is 
that we only need the uniform C°-regularity. We postpone this proof in Appendix 
for self-containedness. 

3.5. Proof of Theorem 3.1: The Regularized Nonlinear Free Boundary 
Problem. We now show the existence of a solution to the regularized free boundary 
problem. 

Lemma 3.10. For each e G (0,eo) '^iT'd S G (0, (^o) with some Eq > and Sq > 0, 
there exists a solution (p^'*^, r^'*^) G C^^^^^ (0^'"^) x C^+"i ([-f , 6'i)) to the regularized 
free boundary problem (j3.2p - (j3.7p . 

Proof. For the notational simplicity, we suppress the (e, (5)-dependence in the proof. 

For each r{9) G IC'^^ C C^+^i ([0^, ^i]), using the solution p of the nonlinear fixed 
boundary problem ()3.2p and ()3.6p - ()3.7p given by Lemma 13.51 we define the map J 
on /C, f = Jr, as in (13. ip : 

(3.36) r{9) = ri + / g{r{s), s, p{r{s), s))ds. 

There are two cases for the approximate shock position f{9): 

Case 1: r{9w) > c(po) + ^- We check that J maps /C into itself. It is easy to 
check that f{9) G C^+^/y {[9^,9i]) D C^+^{[9^,9i)), from the definition of r and by 
using Lemmas l3.3H3.4l property (Ki) follows from (I3.36p . By the definition of g 
and p(-P2) = P, f'{9) = holds, which implies property {K2). Then it suffices to 
show that property (^^4) holds, since the upper and lower bounds of p. Lemma 13.81 
and {K4) imply (-f^s). From the expression of g{r{9),9, p{r{9),9)) and the upper 
and lower bounds of p, we have 

Case 2: r{9w) < c{po) + 5. Since r'{9) > for G {9y,,9i) and ri = c(pi) > 
c(po) + 6, there exists a unique 9a G {9w,9i) such that f{9a) = c(po) + 5. Now, 
choosing r to be determined later such that f{9a + r) < c(po) + 26 and letting 
xi = 9a + T — 9w, modify the approximate shock position on 0^ < ^ < + r by 
defining 

m = c(po) + 6 + Ai9- 9^f + B{9 - 9^)^ 
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with 

A = - —^{na-bxi), B = - -—(6x1 -3a), 

(n — 3)xf [n — 3)x" 

where a = f{9a + r) — c{pq) — 5 and b = r'{9a + t). 
Choose r small enough such that 

bxi — 3a > 0, 

and then n sufficiently large such that 

na — bxi > 0, 

where n depends on 6, but independent of the iteration. Next, we choose n and r 
precisely. In fact, it is easy to see that 

\b\ < 4^Vc^(pi)-c2(Po), 



and 



\bxi\ < C{po,pi,9i,ew) := — — r\/c2(pi) - c2(po)(^i - 9^). 

c(Po) 



If 3(5 < bxi, we choose r such that a = 6 and ni = '^'-^"'^^'^^'^"'^ + 1, which depend 
only on 6, po, and pi. 

If 35 > 6x1, letting r small enough, we can obtain new a and b satisfying 3a = bxi, 
where we choose the biggest r smaller than the old one such that 3a = bxi holds. 
Note that bxi > and r{9a) = c{po) + 6. Thus, choosing n2 = 4, we have A > 
and B = 0. 

Let n = max(ni,n2) = n{pQ, pi,9i,9w,S), which is independent of the iteration 
process. Thus, r(9), uniquely determined, is a strictly increasing function on [9^, 9a+ 
t]. Furthermore, we have 

= f'(-|) < f'{9) < f'{9a + r) = r'{9a + r). 

We define 

ior9€[9a + T,9i], 
''^ ' \m for [9^,9a + T]. 

It is easy to show that Jr{9), 9 G [9w,9i], satisfies properties {Ki)-(K4^). 

First, since f{9) G C^+'y^{[9a + T,9i]), f{9) G C°^{[9^,9a + t]), and {Jr)'{9) G 
C([^«M^i])) we have 

Jr(^)GCi+^^([0^,0i]). 
Next, for 9 G [9^, 9a + r], r'{9) = 'iA{9 - 9^f + nB{9 - 9^Y^^. Then 

r'(02) - r'{9^) 

= 2>A{92 - 9^f - 2.A{9^ - 9^f + nB{92 - 9J^-^ - nB{9^ - 9J^-^ 

n-2 

= 3Ai92 - 9s){92 + 93- 29^) + nB{92 - 93){Y, ^^2(^2 - e«,)"-'-^(^3 - OJ^). 

j=0 
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Using the fact that O2 — Oy^^O-^ — 9^ < xi, and A,B we obtain 

\r{02)' -r{e^)'\ < 1^2-^31" (6^x2-° + C7(n)Sx^-i-°) 

< C(n)(o2;^^-" + 6x7")|02-^3r- 

Notice that f' = |\/ — c^, r G C"'^"'"'^^, and Oa + r is uniformly away from 61, which 
means p £ C^"*"^. We obtain 

which imphes 

\r{92) - r'{e3)\ < C{po, pi,e,S)\92 - if « < ^• 

Thus 

||Jr||ci+-([_|,0i]) < C{pi,p2,e,5) 

if a < min{'yv, |}; which satisfies {Ki)-{K4). 
Thus, we define a map 

J : IC''^ IC''^ 

by 

r = Jr. 

Obviously, IC^'^ is a convex and closed subset of the Banach space C"^, and J is 
compact, if ai < min{^y, ^}. In order to use the Schauder fixed point theorem, we 
need to prove that J is continuous on /C^'*^. 

Assume that r-m, t G /C*^'*^ for m = 1, 2, • • • , — ?^ as jtt, — )• co, and pm solves the 
fixed boundary problem for for each fixed m. Then, by the standard argument 
as in [5], pm — >• p, which solves the problem for r. Therefore, we have 

9{rm{e),e, pm{rm{e),e)) {r{e),e, p{r{e),e)) m ^ oo, 

which implies Jr^ —)• Jr as m —t- 00 at the point where (j3.3p holds for both and 
r. Then Jr^ —)• Jr as m —)■ 00, if Jr belongs to Case 1. For Case 2, due to the 
construction, we divide it into three subcases: 

35 < hxi; 35 > bxi; 35 = 6x1, 

where h = f'{9a + t), xi = Oa ~\~ T and f{6(i ~l~ t) — c(/?o) ~l~ 25 depend only on 

r and 5. For any case, it is easy to deduce that 

{Tm,6a,m) ^ {T,9a), [Am, Bm) ^ {A, B) m ^ OO. 

Then Jr^ — Jr, with the fact that 

Jrm = c{po) + 5 + Ara{9 - O^f + - O^Y 

for < 6a^m + Tmi where n, 9^, and pQ are universal constants. 

Then, for any fixed e, 5 > 0, we obtain the existence of a solution {p^'^ ,r'^'^) to the 
free boundary problem by the standard fixed point argument. Moreover, we have 
^£,5 g ^i+aQgi^^ ^ ^ rpj^-g completes the proof. □ 
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3.6. Proof ofTheorem l3.lt Completion. We note that Lemma l3.10] implies that 
there exists a solution (p^'"^, r^'^) such that r^'^ G JC^'^ . From Lemma 13.61 and the 
interior Schauder estimate, we note that ||/9^'''||^2,q < C, and p^'^ satisfies property 



loc 



(j3.9p . By Lemma |3.7| we have c'^{p^'^) > r^. This completes the proof of Theorem 

EH 

4. Proof of Theorem 2.1: Existence of Solutions 

In this section, we study the limiting solution, as the elliptic regularization pa- 
rameter e and the oblique derivative boundary regularization parameter 6 tend to 
0. We start with the regularized solutions of problem (|3.2p and (j3.5p - (j3.7|) . whose 
existence is guaranteed by Theorem l3.1[ Denote by p^'^ a sequence of the regularized 
solutions of the boundary value problem. 

We first construct a uniform lower barrier to obtain the uniform ellipticity in any 
compact domain contained by r2\(rsonic U rgQ^j^,) for the solutions of the regularized 
problems. 

Lemma 4.1. There exists a positive function ip, independent of e and 6, such that 
(f^O as dist((e, r?), Psonic U T^^^jJ ^ 0, 

and 

c\p''') - ie + v')>^ in n\(rsonic u r^„„i,). 

Proof. For the notational simplicity, throughout the proof, we write p = p'''^. 

1. We claim that c^(p) — (P{p,po) is increasing as a function of p for p > p > pQ. 
In fact, we set 

g{p;po) = c^{p) - c^(p,Po) for p G [p,pi]. 

Then 

P"{p)ip - Po)"^ - P'{p){p - Po) +P{p)-p{po) 



9'{p;po) 



{P-Po? 

In order to show g'{p; po) > for p > p > po, we define 

Hp;Po) =p"{p){p- Pof -P'{p){p- Po) +p{p) -p{po), 

which is the denominator of f'{p)- Then 
(4.1) 

h'{p-.Po) = P"'{P){P-P0?+P"{P){P-P0) 

= {P- Po)((7 - 1)(7 - 2)p^-=^(p - Po) + (7 - 1)/'^-') 

= (7-l)(/0-/Oo)p^^^((7-l)p-(7-2)po) >0 forp>pi>po, 

where we have used the fact that (?{p) = p'{p) = p"^'^. Thus, h{p; po) > since 
h{po) = 0, which implies that g{p; po) > for p > p > po. 

2. Since p > p in fl, we have c^(p) > c^(p) > c^(p, po)- Let 6i be a positive 
constant: 

di := c^(p) - c^(p, Po) < c^(p) - c^(p, Po). 
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Hence, p satisfies c^{p) - (^^ + t?^) > 5i for r = ^JW^^ < in 17. It follows 
from p.2ip in Lemma 13.61 that there exists a constant cq > independent of e such 
that 5i > eo > 0. From now on, we only consider the set: 

n = {P€ n\rsonic : r > ^}. 

For < < 1 and Xq = (^o, %) G ^, let 

CiX) = 1 - + for BniXo) n p.... = 0. 

We define 

where Sq and r are two positive constants to be specified later, and < ^ 1 
independent of e. 

3. Multiplying the equation Qp + eAp = by (7 — l)p'^^'^ and letting c^(p) =: u, 
we have 

22 22 
(4.2) Lu = ^ aijDijU + ^ ai{Diu)'^ + 03 ^ DiuDju + ^ fejAii = 0, 

where 

an = c^(/>) - + e, 022 = c^(/o) - J?^ + e, 012 = 021 = -^??, 

"1 = 1- ^ ^nl^ - + e), a2 = 1 - . ^."^L J ^'(^) - + e)' 

(7-l)c^(/9) (7-l)c^(p) 

(7 - l)c^(p) 

4. Define v = u — (.^^ + 77^) — Using Ln = 0, we have 

= Lu 

2 2 

i,j=l i=l 

2 2 

+03 ^ Alt' + + 7?^ + 95)^j(t^ + + + + 5^ biD,{v + + r/2 + (^) 

j,j=l i=l 

= Liv + 
where 

2 2 

i,j=l i=l 

+4:aiS,D^v + 2aiD^vD^ip + Aa2iDrjV + 2a2Dy^vDy^(p 
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and 

2 2 

^29^ = ^ aijDijip + ^ [ai\Diip\'^ + 6iA<^) + 401^1)599 + Aa2'nDr^ip 

i,j=l i=l 

+2(aii + a22 + 2(^^011 + r/^asa) + , ^^^7 2) " + ^'))- 

(7-l)c^(p) 

We now evaluate L2(p: 
2 

= 5^a,,T(r-iA^C + C"'(T-l)|ACl') 

2 
i=l 

+2 (an + a22 + 2(^^011 + r^\22) + ^^^^^ g'r?^ - 4(e' + r]^)) . 

Noting that, for 

c^(p, Po) 

^-^0 2(7-2)c^(pi) |' 
I (7-l)P0 I 



we have 



2(an + a22 + 2(e^an + ^^a22) + " 2(e' + r/^)) 

P 

> c\po) > 0. 



> 2{e{l - 2(7 - 2)^-t^) + + + 2(0^ - - ry^)) 



Also, since po < p < p < pi, and DjC and Da^, i = 1,2, are bounded in (l, we 
find (52 small enough such that 

L2<f >^c^{po) > for5<S2. 

5. In addition, if i?/j(Xo) n Tshock 7^ 0> we have 

= (7 - l)p^~^Np = Nu = N{v + f + v'^ + ip) 
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and, in O = {c^ — - if < f} C i?j^,(Xo), we have 
N{e + V^) = 2e/3i + 2?7/32 



< 



{-VC + r]){C + m'){c\c' - c2) - 6C{c' + 3c2)) 



Thus, in O := {c^ - - rf < if] C. Br{Xq), by choice of S < 63(61) sufficiently 
small, we obtain 

N{e + r]^ + ip) = 2^/3i + 2r?/32 + 5r(/3i) 

< - V'm + Vv'){c\c' - c?) - 6C{c' + 3c2)) + 6TC~'m\o 

< 0. 

6. Finally, we cover Fq H ^2 with Br centered at Tq so that Dip ■ v = 0. Hence, in 
O C Bji^Xq) with ^0 < min{(5i, 62,63} which depends on p and p2, but is independent 
of e and r < 2, we have 

Liv = —L2(p < 0, 

Nv = -iV(C^ + + V?) > 0, 

Dvv = 0. 

Hence, we apply the weak maximum principle on O C Bpi[Xq) to obtain 

(4.3) inf (c2(p) -f-rf-^)> i^Uc^p) - f - r^^ - = 0- 

O aO 

7. By piecing together these Bsn (Xq), Xq G r2\rsomci we obtain a local uniform 
lower barrier of (?{p^) — (^^ + ^^)- That is, 

c\p) -e-V^>V = 60C in Bs (Xo) n W„ 

where 60 and r are independent of e (though they may depend on R). Moreover, 6q 
tends to as dist((^, r/), Fsonic) — ^ 0, so does y?. This completes the proof. □ 

The proof of Lemma 14.11 also implies that we can obtain the uniform ellipticity of 
(|3.2p which is independent of e in Bsn (Xq) n O^. 

The uniform lower bound of (?—^^—rf' independent of e implies that the governing 
equation (j3.2p is locally uniformly elliptic, independent of e and (5, which allows us 
to apply the standard local compactness arguments to obtain the limit p locally in 
the interior of the domain. 

We first consider the behavior of shock position r^''', as e and 6 tend to 0. We 
divide the shock position into three cases: 



52 ) 



Case 1: c{po) < r{9) < c(pi,po) for all 6 e [6*^, 6*1) and r'{6) = r^J■ 
Case 2: r(6',^) = c(po) and c{pq) < r{e) < c(pi,po), r'{e) = r^J^^^ for all 
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Case 3: There exists 6a G (6*10, 6*1) such that r{6) = c(po) for G [^^.,^0], r{9) > 
c{po), and r'{e) = for 9 G 

Lemma 4.2. T/iere exis^ functions r{e) G C^([0^, 6I1]) and p G Cf+"(ri) n C(n), 
satisfying one of the three cases stated above, such that 

and r) is a solution of the free boundary problem (|3.2|) - (|3.7p . 
Proof. For e, (5 > 0, it follows from Lemma 13.101 that 



r 



where C is independent of e and 6. Thus, by the Ascoli-Arzela theorem, there exists 
a subsequence converging uniformly to a function r{9) in C°([0^,0i]) as e,(5 — )• 
for any a < 1. By the local ellipticity (cf. Lemma l4.ip and the standard interior 
Schauder estimate, there exists a function p G Cf^"^ such that p^'^ — t- p in any 
compact subset contained by il\(rsonic U Lshock); satisfying Qp = in f]. 

For (r(0o),^o) G Tsock with riOo) > c(/?o), there exists a neighborhood of ^0 and 
(5* > independent of e and 6 such that r^''' > c(/9o) + 5* for e and 5 small enough. 
It follows from c{p^'^) > r*^'^ > c(po) + (5* that 

Thus, we obtain the uniform ellipticity locally, as well as uniform negativity oi P -v 
locally. Hence, we can pass the limit to obtain p G C^"^" and 

Mp = on Fshock near (r(6'o), 6*0) 



such that r' [9) = l-^/r^ — c^. 

Then the remainder is to show the case that {r{9o), ^0) G Fshock and r(0o) = c(po)- 
First, it follows from Lemma 13.91 that 

c{po) < r''\e) < c{p''\r'''{9),9)) < c{p''\r''\9o),9o)) 

for 9 G [9^ , ^o] ) and 

c{po) < c{p''\r''\9o),9o),Po) < r''\9o). 

Thus, 

p'''{r''\9o),9o)^ po. 
Therefore, r{9) = c{po) for 9 G [^^,^0]- 

Next we prove the continuity of solutions up to the boundary where r(9) = c(po)- 
First, we prove that r G C^. Still from Lemma 13.91 we easily obtain that 

p{r{9),9) — 7- po if ^ — 7- ^0 from the right. 

On the other hand, 



r'{9) = r{9)\r-^^^^^ foi 9 > 9o, 
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which impHes that r'{0) — )■ as — )■ from the right hand side, and it holds from 
the left hand side obviously. If we define r'(0o) = 0, then r G C^. 
Note that the equation for u = c^{p) is 

(4.4) Q{u) := {c^ - + e)urr + — ^ugg + y '-{urf 



J.2 _ 

1 , ,2 c2-2r2 + e 

+ (^-l),2 M + 1 



= 0. 

We prove the most complicated case ^0 = first, and the other cases will be 
discussed later. 

We construct a family of barrier functions {^'r} with parameter r. For any m > 0, 
there exists 6i{m) > such that r'{6) < m for \6 — 6a\ < 6i{m). This implies that 

\r{9) - r{9a)\ < m5i{m) for \9 - 6a\ < Si{m), 

where Si^m) — ?• as m — ?• 0. 

Let m < 1 and m6{m) = § (t will be specified later). We have 

Po < p{r{e),e) < p{r{ea + di{m)),ea + 5ii7n)) < {c,,r\r{ea + 6iim)))<po + ^. 
For £, 5 small enough, we obtain 

Po<p''Hr'^HO),e)<po + Cm 

and 

< r^^\9) - c{po) < T for \e - da\ < Si{m), 
where C depends only on 7 and pQ. We define 

= = c\po + Cm) + A{c{po) +T-rr + B{9 - Oaf 

in 

Q^'^ = {(r,0) : \r-c{po)\<52^e-ea\<6l{m)}f^n'^\ 

where 82 > t will be chosen later. 
Choose 

R c^(Pi) -c^(Po) . . c2(/ji) -c2(/3o) 

ii = 727—^ , A = Ai- 



5i{m) ' " 5^ 

*r > c2(p^'^) on ag" 



Since po < P^'^ < Pi and p"^'^ < Po + C*"! on F^^'^ , n dQ^'^ , we have 



Thus, we have 

Q(r/) = Aa(a-l)(c2(p^'^)-r2)(c(po)-r + T)"-' + ^ 
(4.5) +(1 - j^ic^ - r^))A^a\c{po) - r + r)^""^ 

"(7-l)r2 



+ (^_lVi^ (^-^a)^ -{c{po)-r + T) . 
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Consider ^ in Q^'^ n {{r,6) : c'^{p''^) - "^p^ > 0}. Since 

c(po + Cm) > c{p^^\r^^\e),e)) > r^'^(0) > r, 

we have 

c^{p''^)-r^ + e > c^{p''^)-^f + c\po + Cm)-r^ + A{c{po) + T-r)'' 
> ^(c(po) + T-r)". 
For a < 1, (|4.5p implies 

< A'a[{2-j:^{c^-r^))a-l){c{p,)-r + ry 



, 2q-2 



-^^Aa{cip,) - r + ry-' + j^B^iO - + 2B^. 

Moreover, let 

A>A, B = ^^. 

If a < 2+c(po pi 7) ^^'^ ^2 + T small enough, we have 

Q^'') < C{po,pi){i2 + C(po,Pi,7))a - + r - r)'"^' + ^^j^- 

0f(m) 

Then there exists a constant A2('52) ""t-, PO) Pi) such that 

Q(^e.'5) < for ^ > A2. 
In fact, if r < c(po)) we choose ^2 = ■v/m'^i("i) to obtain 

c(po) — r + r < 2i/m(^i(m), 

and let 

^(1) ^ C(/?o,pi,a)m^ ^ ^ ^(^0,^1) 
^ <5f(m) ' ^ mt<5f(m)' 

If r > c(/9o), 

c(po) + r - r < T, 

and we let 

(2) ^ C(po,pi,a)?n^^" 
Set = max{Ai, Then 

Passing to the limits 5, e — )• 0, we obtain 

p < ^T- in the domain Q{m, ^/m6l{m)) := n5^£>oQ^'^. 

With these barrier functions, we can show that p is continuous at {r{0a),0a)- In 
fact, for every ei > 0, there exists m > such that 

c'^ipo + Cm)-c^{po) < ^. 
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For this m, we can choose A, B, and r such that 

c'(p) < < y + c2(pi) + A{c{po) - r + r)" + 3(6 - Oaf. 
Choose the neighborhood of {r{9a), 9a) small enough so that 

A{c{pq) - r + r)" < A{2tY < C(po, Pi, a)m?. 
Then, choosing m small again, we have 

c\p)<^ + c\po) + B{e-6af. 

Finally, we choose a small neighborhood such that 

c\po)<c\p)<ei + c\po). 

Thus, we obtain our claim that p is continuous at (r{9a),0a), that is, the results 
hold for this case. 

As for the case £ [9^, 9a), we can choose arbitrary r > 0, which is independent 
of the neighborhood of 9. This fact makes the similar proof of this case much easier 
for all sufficiently small e and 6, and we omit the details here. □ 

Next, we discuss the wave strength at the sonic circle r = c(po) and conclude that 
Case 3 in Lemma 14.21 can not actually occur. 

Lemma 4.3. Let r{9) be monotone increasing in 9 on Fshock o^nd p > po in the 
subsonic region. Then r{9) > c(po) for 9^ < 9 < 9i. 

Proof. We divide the proof into five steps. 

1. We show our claim by contradiction. More precisely, if there exists 9 such that 
r{9) = c{po) := cq. Then, using the monotonicity of r{9), 

r{9) = Co for 9^<9 <9. 

2. For 9o G [9u,,9], we define 

wi = cl + Aiico - r)5 - Si(co - r)^' + Di{9 - 9o)\ 

where Ai,Bi,Di > and ^ < /3i < 1, all of which will be specified later to proof 
that p G C2 near this boundary point. 

Using (j4.4p with the coefficient of Urr replaced by u — r^, we have 

Q{wi) = (-(/32-i)AiSi(co-rf^-i+Oi) 

+ (/3i(2/3i-l)i?2(co-r)2ft-2 + o \ 

(4 6) ^ ^ 

-j^Alic^-r^)ico-r)-^ 

+ {-^{co-r)~l{9-9o)^ + 0,), 
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where 

Oi = -4ico(co-r)-5 -2co/3i(/3i-l)5i(co-r)''i-i + ^(co-r)^ 

+/3i(/3i - l)i?i(co - r)/5i + ^Di + ^(co - r)"^ - /3ii3ir(co - rf'-^ 
+ ^^\t%fH c' - r'){co - rf^-l - ^(c^ - r'){co - r)^ 
+^(c2-r2)(co-r)/5-i, 

O4 = /3l(/3l-l)i3lI)l(co-r)ft-2(^_g^)2 + _4Di_^(^_^^)2_ 

Notice that there exists < a < ^ such that — < (cq — r)" for cq — r > 
smaU. Thus, 

^^^Alic' - r'){co - r)-i| < C{po, pi)Al{co - ry-K 

We can choose a proper constant a such that 

3 1 
/3i--<Q-l, i.e., «>/?!--. 

On one hand, let cq — r > be small enough so that 

{^1 - \)A^B^{c, - rf^-'i > 3C(po,Pi)^?(co - r)"-\ 

which implies 

(4.7) B, > ^^i^°'f^ ^i(co-rr-^-+^ :=AMPo,Pi,M{co-rr-f^^^l 
Pi 4 

On the other hand, if cq — r > is sufficiently small, we have 



which implies 



(4.8) A > ^^^^f#^(co - r)^-i := C{P,)B,{co - r)^^~l 

Pi ~ 4 

Moreover, we have 

C{po,Pi,l3i){co - rr-^^+'2 < C(/3i)(co - r)^"^^ 

when r G [f, cq] and f is close to cq. 

Choose proper constants >li and Di such that 

u;i > eg + ^Ai{co - r)i + i?i(^ - ^o)^ > 

at the boundary of a relatively neighborhood A'^i of {co,Oq) to $7. Choose i?i suffi- 
ciently small such that 

Q{wi) < 
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and 

C(po,Pi,/3i)(co - rr-^+-^ < 1^ < mm{c(/3), ^^^^-l^}(co - rf^-^^ in N^. 

This implies that (02]) and (jMI) hold. 
Obviously, we have 

drrWi < in A''i 

if (HZ!) and (USD hold. If 5i = {{r,e) G iVi : > u;!} / 0, we have Q{wi) < 
Q{wi) < in ^i. Thus, 

< Qu- Q{wi). 

Using the maximum principle, 

which contradicts with (? > wi. Thus 

(? < wi in ^1- 

3. We define 

W2 = cl + ^2(co - r)5 + B2(co - rf^ - D2{9 - Oof, 

where A2, B2, D2 > and ^ < /32 < 1, all of which will be specified later to prove 

that C2 is optimal. Through a simple algebraic calculation, we have 
(4.9) 

QM = m -l)A2B2{co-rf'-UOi) + {(32(2(32 -l)Biico-rf ^'-' + 02) 



+ ( - (32i(32 - l)B2D2{co - rf^-\6 - Oof + O 



+ {(32B2r{co - rf^-^ + O4) + i^2l)2(co - r)^2(0 - O^f 
+2P2{P2 - l)B2Co{co - rf^~\ 

where 



2 



r 



2 ^ ' 2 ' ' 2r 



3 



2 



(7 - l)c^ 
(7-2)/3ii?| 

(7 - l)c^ 



- -(^(^-^0)^ 
O4 = -ME(<;^_r^)(e„-,.) 
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Let D2 be large enough such that > W2 for some 6 = 6a, Ob- We choose r < cq 
such that 

c2 > d + 2A2{co - f)^^ - D2{e - 60)^ 

> CO + ^2 (CO - r)H B2{co - ff^ -D2{9- O^f . 

The second inequahty holds provided that ^ < (co — f)2~^2_ Choosing /32 > |, we 
have 

\hB2r{co - rf'-^ + 2W2 - l)B2Co{co - rf^'^ < for ^ < r < cq, 

and 

Q{W2) > 0. 

Then, if ^2 = {(r, 9) e Ni : < W2} / 0, we have 

Q{w2) > Q{w2) > in 52. 

Thus, 

Qu - Q{w2) < 0. 

Using the maximum principle, c < W2, which contradicts with < wi. Thus 

> W2 in A''2. 

4. We now show that 

c^>cl + AsicQ - r)i + ^3(00 - r)^3 ^3 

in a relative neighborhood of {ro,9o), where A3 and B3 are positive constants to be 
specified later, so that the C 2 -regularity is optimal. 
Since (P > W2, we can choose 9a and 9h such that 

> eg + A2{9a,9a){co - r)^ + B2{co - rf^ for C iV2. 
Thus, there exist positive constants j43, i?3, and ^-^ such that 

U'3 < c^. 

It is easy to see that 

QM = ((/3i-^)^3i?3(co-r)'^^^i+Oi) + (/33(2/33-l)i?|(co- 0^^^-2 + 02) 
+ (/33i?3r(co - rf'-^ + O4) + - l)i?3Co(co - r)/^3-\ 
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where 



4(7-l)c2' 
H — ^(co-?^) "^(c -r)(co-r 

(7-2)/33A3i?3, 2 2 



(7-l)c2 



02 = /33(/33 - 1)^^3(00 - r)fe - ^^^^ (c^ - r^)(co - r)^fe-2^ 

03 = -^(c2-r2)(co -0/3-1. 

r 

Similarly, we can show that (? > w-^ in N-^. 
Thus, we have 

^A3(co - r)^ < - eg < 2^i(co - r)i in A^i n iV3- 

This implies 

a(co - r)5 < 77 := /9 - po < ^(co - r)5 in A''i fl A''3 
for some constants a and ^, so the optimal regularity of p is C2 near the sonic circle. 
5. We introduce the coordinates: 

= (co - r, 6* - 6^) 

and set 

2 2 

V = C — Cq. 

Thus, rewriting the equation for (? in the divergence form, we have 

(4.10) Qv = {aii{v + 2cqx - x'^)vx)^ + hv^ + ia22Vy)y = 0, 

2(2-7) 2 2 

where an = 022 = (Tpjp-, and 61 = 

Scaling v in Ni D by defining 

(4.11) uiS,T) = ^viS-'i,yo + S-'^T), 

12 14 

for {S~~5 ,yQ + S^~T) G Ni n N3. Moreover, u satisfies the following governing 

equation: 

(4-12) 

Qu = iauus)s + iai2UT)s + {a2ius)T + {a22UT)T + (^2'w)t + hus + huT + d2U = 0, 



SHOCK DIFFRACTION BY CONVEX CORNERED WEDGES 



39 



where 



^11 = aii{S^u + 2cQS 5—S 5), 

147^ 7 6 18 

ai2 = -=-^au{S^u + 2coS~5 - S~~), 

14T 7 6 ri~^\ 

0-21 = -TTrOiiK!^'''U + Icob s-i 5), 

144 189r2 ^7 ^ ^„i8, 
022 = — «22+ ^^g,2 aii{b^u + 2cQS 5-S 5), 

h = -^^au{S^u + 2coS~^ - S~^) = b22TS"'^, 

S'^U(2-J)S~^ 2(3-27) 7 ^ ^_6 ^ 

5(7 - 1)^ 

12 

aii(4co-S'"^) 1261 _ n 

n IT = cii-^ 5 , 

55- 55- 

168r I 6 M 189ran IGS^iT . 2 

7 16 
~ S5u{2-'y)S~^ 2(3-27) 7 6 18. 

^ = — ^^7 T\2 — c 55-u + 2co5 5-S 5 

25(7 — 1)^ 



18 , 



Cl 



13aii(2co - 12aiiS'in(co + f5-t) 1261 7 ^_i6 



255*^ 25Sf 25Sf 

From the optimal continuity, 

< a < S^u < A, 

we have 

< C'^ < Ai, A2,622,Cll,C22,C^l < C 

if and T are sufficiently small. Here Ai and A2 are the eigenvalues of the matrix 
(oij)2x2i so the equation is uniformly elliptic for u in the (S, T)-coordinates. 

Let Xq < 5" < Xg with xq small enough. Then, using Theorem 8.20 in [16], we 
have 

clxq < n(xQ , 0) < sup u{S, T) 

3Iq ^ S "-^Q 



< C7 inf ^ ^ n(5, T) < C7u(xo ' , 0) < CAx] , 

3? f 1 *^ 5 3? c 



where C < C(n)'^A+^-^) in [16] is independent of xq, since (A, A) = (Ai,A2), R = 
•^0 * ~ "^0^ — ''^o ^' ^'^'^ ^ max ^ {&2)Ci,C2) V c^} < Cxq . This implies that 

<S<Xn 
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7 



Xq ^° < C, which is a contradiction if xq is sufficiently small. This completes the 
proof. □ 

Next, we consider Fshock in the ?7)-coordinates to obtain finer properties. 

Lemma 4.4. For the free boundary Fshock = {(CiViO) '■ < £, < Ci} determined 
by ([3:2]) - (13771) ■ we have 

r/(e)GC2([e«„6)) 
and r]{^) is strictly convex for G ['?t<;,Ci)- 

Proof. We define 

(4.13) F{^,r])=f+r^^-r\e{^,rj))=0 on Fshock. 

It is easy to check that 

= {2rj - 2rr'a,)|g=5„ = 2r]{^J + 0. 

By the implicit function theorem, there exists r/ = r/(^) such that ()4.13p holds locally 
on Fshock near C = That is, there exists ^ > such that {S,,r]{^)) G Fshock for 

< c < e 

Recah that r?'(0 = f{^,v{0, Pi^^viO))- Then 

V" = h + fW + fpP for ? G (^«,, 0- 
Observe that, if p were constant, the shock would be a straight line. We conclude 

+ fW = 0. 

Therefore, the sign of 77" is determined entirely by the sign of fp and p' . Note that 
p is increasing, p' > 0, and ^ > 0. Moreover, we have 

dl ^ -2^r,c,/e + r?2 - c2 + 2r,\e + - (?) + (g^ + v''){<? - y^) 
d(? c{ir] - c^e + - c2) + r?2 - c2 



,2 



c{^rj — C\J ^2 + 7^2 _ g2) 2 _|_ j^2 _ g2 

If gr/ < 0, it is clear from (|4.14|) that 

— 

If gr/ > 0, from ()4.14p . we have 



dc^ 



C-sJ ^2 + 7^2 _ c2|'^2 _ c2)(7^2 _ ^2^ _|_ jyy^ ^2 _|_ j^2 _ ^2 



These imply that r] = r]{S^) is strictly convex for ^ G [CwnCi)- 1^ 
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Based on the lemma above, we have 

Remark 4.5. Problem (j3.2p - (j3.7p is equivalent to the following free bomidary prob- 
lem in the self-similar coordinates: 

(i) Equation: 

2 2 

(4.15) Lp = Diia^j{C,rj,p)Djp) + Y,h{i,ri)Dip = in 

i,j=l i=l 

with 

ai2 (C, ^, P) = «2i (C, ??, /j) = -ir], bi{^,r]) 62 r/) = ry. 

(ii) The shock equation: 

= f{C,V,p) = ^2 _ c2 with 7?(6) = 7?i, 

with the boundary condition on Fshock: 
2 

(4.16) Np = Y,(^^D,p = on Tshock = {v = ^(0 : < ^ < 6}, 

1=1 

where /3i and /32 are the following functions of (^,r?), p, and r/': 

(4.17) /5i = {f + r^''){-^'^ + r^){c\p) + c\p,po)) 

-2c\p,po){- r?'e(c' + r/^) + iv - viv'f " er?')(c' " f)) 

and 

(4.18) /?2 = v'{e + ^'){v-v'0{c\p) + c\p,po)) 

-2c\p,po){{v'v - e - avT)ic' - v') + ^'^(c' + f))- 

(iii) The remaining boundary conditions: 

(4.19) P = P2 on Tsonic, Pu = on To, /o(P2) = P, 

where v is the outward normal to 0, at Tq. 
It is easy to check that (j4.16p is the oblique derivative boundary condition along 

-Tshock- 

With Lemma 14. 4^ we can show that Case 1 is the only case for the solutions, 
which implies that we can obtain the finer regularity near P2. 



Lemma 4.6. Suppose that (p, r) is the solution to the free boundary problem (j3.2 
(j3.7p . Then the shock does not meet the circle r = ro at the wedge. 



42 GUI-QIANG G. CHEN XUEMEI DENG WEI XIANG 

Proof. We show our claim by contradiction. Otherwise, r{6w) = cq. 

1. First, let ry = rcos(0 — 0^) and consider 

(t> = cl + Ai{co - - Si(co - r^f^ + Ci{e - Oof, 

where 9 G [OwjOw + S], 6 > small enough, Ai,Bi,Ci > and ^ < /3i < 1, all of 

which will be specified later to prove that p £ near this boundary point. We 
now show that (p is an upper barrier of p. 

2. Since < Cq — ry^ = (cq — r^) + r^sin^(6' — 6^) on Fshock from its convexity 
indicated in Lemma 14.41 we have 

0<r^-cl<r^ sm\9 - 9^) < C{9 - 9^f on T^hock 

for some positive constant C . This implies that 

- cl<r'^ - <r'^ - cl< C{9 - 9^f. 

Then 

c^-cl< C{9 - 9^f, 

since and are both functions of p. We can choose Ci > so large that < </> 
on Tshock- Moreover, we have 



wedge ' 



- wedge 



/3iBi(co - r]f-^)rsm{9 - 9^) + 2Ci{9 - 9^) 



where v denotes the unit normal. We can also choose proper Ai and Bi such that 
< in the remaining boundary of a neighborhood of {ry,,9y,) as in Lemma 14.31 



3. Using (j4.4p for Q with the coefficient of u„ replaced by n — r^, we obtain Q. 
That is, 
(4.20) 

:= (-/3ii?i£i(co-7?)/^i-2+Oi) 

-iyli(co - r?)-| (Ci{9 - 9of cos\9 - 9^) + sm^{9 - 9o)) - 402, 
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where 

-/3i(/3i - l)Bi{^ - r2)(co - r?)^i-2 cos" {6 - 6^) 

4 zr 
+ (^ - /3i(/3i - l))AiSi(co - cos2(0 - 0^) 

+/3i(2/3i - l)Bl{cQ - r^ff^-^cos^e - 9^) - (^lA^co - vf'''^ cos^ 
(7 - 2)AjBf , 



and 



(7 - l)c' 



2 



cos 



_^(c2 _ 2r2)(co - 7]r-2 cos{9 - 6^) 
Zr 

(7 - l)c^ 



O2 = - l)BiCi{co - v)^'-\e - OoY cos 



2 „2/ 

■ - ^0) 

-/3i(/3i - l)5ic2(co - r/)^i-2 sin2(0 - 0^) 

7 — 1 (7 ~ Ij'^ 

PiAiB 



2 



7-1 

+ /'^' (CO - r?)-5(0 - 0^) sin(0 - 0^) 
(7- 1)^ 

_ 4/3ii?iCi _ ^ ^^„i _ ^.^ _ 
(7- l)r 

It is easy to see that there exists a small neighborhood N\ of {j^w^^w 

) such that 

g((/.) < in A^i. 

Obviously, we have drr4' < in A'"i. 

4. Now we show that c^{p) < (p. If = {{r,6) £ Ni : > / 0, we have 

Qiipi) < Qii^i) < in Si. 
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Thus, Q{u) — Q{il^i) > 0. Using the maximum principle, we have u < ipi. This 
contradicts with > V'l- Thus, < V'l in ^i, i-e., 

< c2 - < Ai(co - r/)5 + Ci{9 - e^f. 

5. Now we find a lower barrier for p. 

As the proof of Lemma 14. 3| we can show that there exist a neighborhood N2 of 
{fw^(^w) and a positive constant A2 such that 

- eg > ^2(co - r)5 in iV2 n {(r, 0) : r < cq}. 

The only new here is the boundary r = cq, which is obvious. This implies that 

a(co - r)^ < ?; := p - po < -4(co - r)5 in iVi n iV2 n F, 

where V is an upward sector containing the wedge, with the vertex at P2 and the 
angle smaller than ^, for some constants a and A depending on V , so the optimal 

regularity along the wedge is C2 near the sonic circle. 

6. We introduce the coordinates: 

(4.21) X = Co - r, y = e - Ow, V = - cl- 
Thus, rewriting the equation for (? in the divergence form, we have 

(4.22) Qv = {aii{v + 2cqx - x'^)vx) ^ + biv^ + {a22Vy)y = 0, 

2(2-7) 2 2 

where an = " ^'^'l , 022 = t^^, and 61 = 
Scale f in iVi n A/'s n y by defining 

(4.23) u{S,T) = ^v{S-f,S~^T) 

S 5 

for (S^^jS'^^T) G A''i n n y. Moreover, m satisfies the following governing 
equation: 

(4.24) Qu = 0, 
where 

(4.25) Qu := {aiius)s+iai2UT)s+ia2ius)T+ia22UT)T+{b2u)T+ciUs+C2UT+d2U, 
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with 



^11 = aii(5'5 n + 200^ ^—S s), 

14T 7 6 18 

021 = --^aii{S^u + 2coS 5-5 5), 

ob 

144 189r2 ^7 ^ ^_i8, 

022 = -^a22 + -^^^an{S^u + 2cob s-b s ), 

14T 7 _6 _18,~ _2 

b2 = ^g^aii(5'5u + 2co5 5-5 5 ) = 622TS' , 



Cl 



55tl(2-7)S'-- 2(3-^7) ,^7 „ ^ 6 ^18^ 

5(^_1)2 g (5iu + 2co5-5 -S"^) 

aii(4co - S'^) 1261 , 
n IT - Clio 5 , 

168T 1 6 18 189ran 168hT . 2^ 

~ slu(2--f)S-f ^(ft±d,7 ^ ^_i8^ 
= — 25(7 -1)2 — c h-i) (55u + 2co5 5-S 5) 

13aii(2co — S"^^) 12aiiS'5 n(5co + 65" 5 ) 126i „_16 
16 23 16" = "i*^ ^ • 

255*- 1255- 255*- 

7 

From the optimal continuity, 0<a<S5u<A. Then we have 

< C"^ < Ai, A2,622,Cll,C22,C^l < C 

if and T are sufficiently smah. Here Ai and A2 are the eigenvalues of the matrix 
{o-ij)2x2, which implies that (j4.24p is a uniform elliptic equation for u in the {S,T)- 
coordinates. 

Let Xq < S < Xq* with xq small enough. Using Theorem 8.20 in [TB], we have 

7 

axQ <u{xq^,0) < sup u{S,T) 

-5/4 

< C inf u{S,T) KCuix'f.^O) <CAxL 

-5/4 
[ 

where C < C{n)^'^^^ in jl6j is independent of xq here, since 

(A,A) = (Ai,A2), R = x~~' -x^^ <x~\ 
11= max {b2,ci,C2,Vd} < Cxq. 

''^ iS* "-^Q 
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This imphes that 



< C7, 



which is a contradiction if xq is small. This completes the proof. □ 

Finally, we establish the Lipschitz continuity for the solution near the degenerate 
sonic boundary. 

Lemma 4.7. The solution p to the free boundary problem ()4.15p - ()4.19p is Lipschitz 
continuous up to the boundary T sonic- 

Proof. On one hand, since p < pi in il., we have 

c\p)-e-v'<cHpi)-e-^f. 

On the other hand, it follows from Lemma 3.1 that 

c\p) - f - ri^ > f + v"^ - c\pi) inn. 
Let r| = (?{p2)- We have 

\c\p)-c\p2)\ < \c\p)-e-ri^\ + \c\P2)-e-n'\ 



< 2\c\p2)-e-v 



,2 1 



which implies that p is Lipschiz continuous up to the degenerate boundary rgonic- D 



Proof of the Existence Part of Theorem 12. ll The above seven lemmas, i.e., 
Lemma l4.H -Lemma 14.71 show that there exists a solution 

which satisfies ()2.1ip - (|2.15p . This completes the proof of the existence part of 
Theorem O □ 



5. Proof of Theorem 2.1: Optimal regularity near the sonic boundary 

In this section, we prove that the Lipschitz continuity is the optimal regularity for 
p across the sonic boundary Lsonio as well as at the intersection point Pi between 
Tsonic and Fshock- In §4, we have shown that the solution p to the free boundary 
problem ()4.15p - ()4.19p is Lipschitz continuous in 0, up to the degenerate boundary 
rsonic- Now we employ the approach introduced in Bae-Chen-Feldman [2] to analyze 
the finer behavior of p near the sonic circle r = ri := c{pi). 

For e G (0, we denote by 

:= n n {{r,e) : < ci - r < e}, 

the e— neighborhood of the sonic circle Fgonic within Q. 
In r^e, we introduce the coordinates: 

(5.1) {x,y) = {ci-r,e-ei). 
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(5.3) 



CI ■ 



'^{2cix — -0 — x"^), 



One of our main observations is that it is convenient to study the regularity in terms 
of the difference between (pi ) and (p) : 

(5.2) iP:=c\p,)-c\p), 

since ifj and p have the same regularity in il^. 
It follows from (|1.13p that satisfies 

{2ciX - V + Oi)^^^ + (Ci + 02)V'a; " (1 + 03)0^ 

+(1 + 0,)^^;yy - ( + 05)^2 ^ 
in the (x, y)— coordinates, where 

'01(3;, V') = -x^, 
02(2;, 0) = -Sx + 

(5.4) i03{x,ip) 

0^{x,^) 
Moreover, ip satisfies 

(5.5) > in 
and the following Dirichlet boundary condition: 

(5.6) -0 = on n {x = 0}, 
where 

(5.7) := {(x, y) : X G (0, r), |y| < C M^, 

with R = 9w—di, since we can extend ^(x, y) from 0,^, by defining ?/;(x, y) = ^(x, — ?/) 
for {x,y) S ilj, and extend the domain il^ with respect to y. Thus, without further 
comment, we study the behavior of ip in Q^^- 

It is easy to see that the terms Oi{x,y), i = 1, • • • ,5, are continuously differen- 
tiable and 



(ci-x)'^ 

1 

(ci-x) 



A7 



h 



1 



(5.8) 
(5.9) 



|Oi(x,y)| |Ofc(x,y)| 



\DOi(x,y)\ 



< N 



X 



for k = 2, - ■ ■ ,5, 
, \DOk{x,y)\<N for A; = 2,--- ,5 



in {x > 0} for some constant depending only on ci and 7. Inequalities ()5.8p and 
()5.9p imply that the terms Oi{x, y), i = 1, • • • , 5, are small. Thus, the main terms 
of (|5.3p form the following equation: 

1 



(5.10) 



(2ci - i>)iJxx + CiV'x + 



yy 



(7 



1)4 







in Q 



r,R- 
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It follows from Lemmas 4.1 and 4.5 that 

(5.11) < V < 2(ci - t?)x, 

where depends only on pi and 7. Then equation (j5.10p is uniformly elliptic in 
every subdomain {x > 6} with 6 > 0. The same is true for (j5.3p in if r is 
sufficiently small. 

Remark 5.1. If r is sufficiently small, depending only on ci and 7, then (|5.8p - (|5.9p . 
and (|5.1ip imply that (|5.3p is uniformly elliptic with respect to ip in n{x > 5} 

for any 5 € (0, §). We will always assume such a choice of r hereafter. 

5.1. First-order lower bound of ip. In order to prove that C^'^ is the optimal 
regularity of ip across the sonic boundary, our idea is to construct a positive sub- 
solution of (|5.3p . (|5.5p . and (|5.6p first, which provides our desired lower bound of 

Lemma 5.2. Let ip he a solution of the Dirichlet problem ()5.3p and ()5.5p - (j5.6p . 

T/ien i/iere exisi r > and ^ > 0, depending only on ci, 7, 0^, and inf 

Qt«n{x>f/2} 

suc/i that, for all r G (0, |], 

(5.12) V'(3;,2/) > inQ^^. 

16 

Proof. In the proof below, without further comment, all the constants depend only 
on the data, i.e., ci, r, 7, 9w, and inf ^p, unless otherwise is stated. 

Qfj,n{x>f/2} 

Fix yo with |yo| < We now prove that 

5 

(5.13) ip{x,yQ) > -fix forxG(0,r). 

8 

Without loss of generality, we may assume that R = 2 and yo = Oi otherwise, we 
set iP{x, y) = iPix, yo + §y) for all {x, y) G Qt^. Then ^{x, y) G C(Qt^) n C2(Q+^) 
satisfies (j5.3p with (j5.8p - ()5.9p and (j5.1ip in Q^2' with some modified constants N , 
"d, and Oj, depending only on the corresponding quantities in the original equation 
and on R. Moreover, 

inf ip = inf ip. 

(9+n{x>f/2} Q+^n{x>r/2} 

Then (j5.13p for tp follows from (j5.13p for 'ip with yo = and R = 2. Thus we keep 
the original notation with yo = and R = 2. That is, it suffices to prove that 

5 

(5.14) Tpi^x^O) > -fix for X G (0, r). 

By the Harnack inequality, we conclude that, for any r G (0, |), there exists 
a = a{r) > 0, depending only on r and the data ci, r, 7; (^wj and inf 

Q+^n{x>r/2} 

such that 

(5.15) ^ > (T on 3/2 {"^ ^ 
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Let r E (0, §) and 

(5.16) < ^fo < mm{^,ci}, 

r 

where r will be chosen later. Define 

iKy + '^?, -i<y<-h 

(5.17) ff(y)= ^i(2/-2y2 + |), -i < y < 1, 

Set w{x,y) = ^xg{y) with g G C^([— 1, 1]). Then, using (|5.16p and (|5.17p . we obtain 
that, for all x G (0,r) and \y\ < 1, 



(5.18) 

Therefore, we have 



t/;(0,y) = 0< V(0,y), 
-(^(r,?/) < |/ur < y), 
_u;(x,±l) = < 'il){x,±l). 



w < ip on dQ'^]^. 

Next, we show that w{x,y) is a strict subsolution Ciw{x,y) > in Q^i, if the 
parameters are appropriately chosen. In fact, 

Ciw{x,y) 

= {cig{y) - g^{y)) 

+x[g"iy) - ^-L^x(5'(y))2 + ^g{y) - ^g\y) + 0,g'{y) - xO^{g' {y)f) . 
On one hand, for 1 — |?/| < Eq with small enough, we can see 



a'iy) - T^-rV2<9^y)f + —9{y) - —a'iy) + 0,g'{y) - xO,{g'{y)f 
(7 — Ijcf X x 

1 



^ 9"iy) - ^^TTTa;(<7'(?/))' - Nx{g{y) + l)g{y) - xNg'{y) + Nx\g'{y)f 

(7 - i)q 

=: h{x,y). 

It is easy to see that /i(x, y) is continuous with respect to x, /i(0, y) = 0, and that 
there exists ri > such that h{x,y) > for r <ri. 
On the other hand, for 1 — |y| > eo; 

£itt;(x,y) 

> x(g"iy) - ^-L^x(5'(y))2 + ^giy) - ^g\y) + 04<7'(y) " xO,ig'iy)f 
+/xeg(ci - 1^). 

Then there exists r2 > such that the above inequality is positive. 
We claim 

sup(it; — il^) < sup {w — iIj) < 0, 
Q+, dQt,i 
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whenever < r < tq := min{ri,r2} and /U G (0, /.to]- Otherwise, there exists a point 
{xo,yo) S Q^i such that 

< {Ciw - Ciip){xo,yo) 

= {2cix + Oi){w - + (ci + 02){w - - (1 + 03){w + i>)^{w - tp)^ 

+ (1 + 04){W - i;)yy - {- ^ + O5) {W + ij)y{w " l/j) y < 0, 

where we have used the fact that Wxx = 0, which is a contradiction. Hence, we 
obtain our claim: 

y) > w{x, y) = xf{y) in Q+p 

In particular, 

5 

ip{x,Q) > -fix for X G [0, r]. 
8 

This implies (j5.13p . Then (j5.12p holds by modifying ft, which is still denoted by fi. 
This completes the proof. □ 

5.2. C^'"-Estimate of ^. If •0 satisfies ([53]), (I53])-(I5S]), and (f5TT]l . it is expected 
that ip is very close to cix, which is a solution of (IS.lOp . More precisely, we now 
prove 



— cix\ < Cx'^'^" for all {x,y) G Qj" 



TR 
8 



for some constant C. 

To prove this, we study the function 

(5.19) W{x,y) := cix - 'ilj{x,y). 
By dSS]), W satisfies 

(5.20) C2W = {cix + W + Oi)Wxx - (ci -O2- 2ciO^)Wx + (1 - Oa)^^' 

+ {l + 0,)Wyy-{j-^^-0,)W'y 

= C1O2 + cfOs in Qfj^, 

(5.21) VF(0,y) = on ^Qt^ n {x = 0}, 
and 

(5.22) - (ci - '&)x < W{x, y) < ax in Qt^. 

Then we establish the following two estimates. 

Proposition 5.3. Let ci, f, R, and "d be the same as in Lemma \5.2[ Then, for any 
a G (0, 1), there exist positive constants r and A, which depend only on N , ci, r, R, 
•d, and a, such that, ifWe C{Qf j^) D C2(Q+^) satisfies (ICTj) - ([5:22]) . then 

(5.23) Wix,y) <Ax^+'' inQ^^. 

^' 4 
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Proof. The main idea of the proof is the same as that in [2], and we only hst the 
major procedure and the difference here. The complete proof is very long, which 
can be found in Appendix. 

First, we prove that there exist ai G (0, ^) and ri > such that, if G C'(Qt^)n 
C2(Qt^) satisfies §M>~<^M), then 

whenever a G (0,ai], r G (0,ri], and fii < min{/i, i}, where fi is the constant 
determined by Lemma l5.2i 

As in [2], we first note that, without loss of generality, we may assume that R = 2 
and uq = 0. Then it suffices to prove that 

W{x, 0) < for X G (0, r) 

for some r G (0,ro) and a G (0, ai), under the assumptions that (|5.2U|) - (|5.22|) hold 
in Qf^ and (jX25]) holds in (5+ 2- 
For any given r G (0,ro), let 

Air = ci{l - m), Bi = ci{l-ni), 
V = Aix^+''{l-y'^) + Bixy^. 

Then we obtain 

w<v onag+i, 

and 

C2V - C2W - v^^{v - W) < inQ+i 
whenever r G (0,ri] and a G (0, ai] so that 

(5.24) (2a-l)(a + l)ciAi < 

and 



(5.25) ,,<„„„nr^j.,^r^j.,,„|. 

Here we have used the fact that 

v,AW -v) = Q(a + - y^)iW - v) < {I + a)aci{l - - 

Noting that 

C2V - C2W - Vxx{v - W) 

= (cix + W + Oi){v - WU - (ci -O2- 2ci03)(,v - W)x 

+ (1 - 03){V + WUV -W,, + {1 + Oi){v - W)yy 

\,2 - 05)iv + WUv - W),, 

(7 - i)q 

which implies 

W <v in g+i. 
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Next, we generalize the result for any a £ (0, 1), which suffices to show that for 
the case a > ai. 

Fix any a G (ai, 1) and set the following comparison function: 



"l^o— oi ^ t) J ' oil a 



2 



Then 

(5.26) W <v on 5Q+i for r G (0, n] 

As before, we can prove that 

C2V - C2W - v^^{v -W) <0. 
Then it is easy to prove that this proposition holds with 

ci(l - m) 



A 



^'^1 j^a—ai 



This completes the proof. □ 

Proposition 5.4. Let ci, f, R, {), and Oj he the same as in Lemma 5.2. Then, for 
any a G (0, 1), there exist positive constants r and B, depending on N , ci, r, R, 
and a, so that, ifW £ C{Qj ^) n C^iQjpt) satisfies ([OO]) ^ (|5:22]l . we have 



(5.27) W{x,y) > -Bx^+'^ in Q 



l+a 

3-R ■ 

4 



Proof. Similar to the proof of Proposition 15.31 it suffices to prove that, with the 
assumption R = 2, 

W{x,0)> - for 2; G (0, r) 

for some r > and a G (0, 02)- For this, we use the comparison function: 

v{x, y) := -Lx^+"(1 - y^) - Kxy^, with Lr" = K = 
It is easy to check that 

W > V on dQ^-^ for r G (0, n]. 

Then we follow the same procedure as in [2], except that C2V > C2W, to find that 
the conditions for the choice of a, r > are inequalities (j5.24p and ()5.25p with 
(/Lii,ri) replaced by (/3,r2), respectively, and with an appropriate constant C. 
We claim that 

min(H^ — v) > min(VF — f ) > 0. 
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Otherwise, there exists a point (xo,yo) £ Qti such that {W — u)(xo,yo) < and 
> {C2W - C2v){xo,yo) 

= (cix + W + Oi){W- vU - (ci -O2- 2ciOs){W - v)^ 
(5.28) +(1 - 03){W + v),{v - W), + (1 + Oi){W - v)yy 

> inQ+, 

which is a contradiction. This completes the proof for the case a < 02- 
For the case a S {a2, 1), we set the comparison function: 

^-(^'2/) -"T^v^^ + )-— ^ 
'2 ' '2 

Then, using the argument as before, we can choose r > appropriately small such 
that 

C2U-. - C2W > 

holds for all G Qti- 

This completes the proof. □ 



Lemma 5.5. Let ip G ^iQtn) ^ ^"^iQt r) ^ solution of the Dirichlet problem 
(OD and (I53D -(15:61). T/ien ^/^ G C'^'°(Q5/2i?/2) " ^ (0, 1) with 

R 

V'x(0,y)=ci, Vy(0,y)=0 for any\y\<—. 

Proof. The proof is quite similar to that in [2] , and the main difference is the scaling 
due to the different equations. Here we just list the main procedures. As in [2], we 
divide the proof into the following four steps. 

Step 1. Let V' be a solution of (j5.3p in Qf^ for f as in Remark 5.1. By Theorem 

12. H we have 1I) G C^'^^iQ^ P112) ^'^^ some r G (0, r/2). Thus, it suffices to show that, 

for any given a G (0, 1), there exists r > so that ip G C^~^'^{Q'^ ^j^) and 

'0x(O,y) = ci, V'j/(0,2/) = for all |y| < |. 

Let y) be defined by (|5.19p . We need to show that, for any given a G (0, 1), 
there exists r > such that 



(i) e 

(ii) L»(0,y) = for all |y| < f . 

Step 2. By definition, W satisfies (|5.20p - (|5.22p . For any given a G (0,1), there 
exists r > so that (fOSl) and (jOTl) hold in Q^^^j^i^ by Propositions lOHOl Fix 
such r > 0. 
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Furthermore, since W satisfies estimate (j5.22p . we can introduce a cut-ofF function 
into the nonhnear term of (j5.20p . That is, we modify the nonlinear term away from 
the values determined by (|5.22|) to make the term bounded in — . Namely, fix 
C e C°°{R) satisfying 

-(ci - |)x < C < (ci + f)x onM, 

(5.29) C{s) = s on (ci - i?,ci), 
C(s) = on (2ci -i?,ci + /3). 

Then, from (3.26) and (3.28), it follows that W satisfies 

x{ci + C(f ) + Oi)W,, - (ci -O2- 2ci03)W,, + (1 - Os)W^ 

(5.30) +(1 + 0,)Wyy - ((:^ - 05)%' 
= C1O2 + cf O3 in Qt^, 

Step 3. For z := (x, y) G i?/2' define 

Rz ■■= {{s,t) : < |, \t-y\< ^}. 

Then 

(5.31) Rz C Q^3R for any z = {x,y) e Qt r- 

'4 2 ' 2 

Fix zo = (xo,yo) G Qt/2,R/2- Rescahng in R^g by defining 

^^^^''^5,^) = ^^(^0 + ^S,yo + ^T) for {S,T) G Qi, 

where Q/^ = {—h,Kf' for /i > 0. Then, as [2], l^(^o) satisfies the following equation 

for (5,r) G Qi: 

(5.32) 

2 2 
^ Aii(VF(^o),5,T)L'2.iy(^«) +^Bi(Z)M^(^o),5,T)AVF(^"^ = in Qi, 

i=l i=l 

and, if r > is sufficiently small, depending only on the data, then ()5.32p is uniformly 
elliptic with elliptic constants depending on ci but independent of zq, and that the 
coefficients Aii{W^^''\ S,T), B,{DW^^°l S,T), and F(^o)(5,r), for p G M^, {S,T) G 
Qi, satisfy 

11^ '''' IIcO{QT) - I^' II^mIIcO(RxQT) - ^' ^ri')llcO(R2xQr) - 

7* 2 

where C depends only on the data and is independent of zq. Then, by Theorem 
(Al) in [7], 

WW^'^^^^^ioTj;) ^ ^(ll^^^"^llco(QT) + l|i^(-)|lcO(gr)) 
< + =:C, 
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where C depends only on the data and a in this case. From (|5.33p . we have 

(5.34) \D^W{xo,yo)\ < Cx^, \DyW{xo,yo)\ < 

for ah (xo,2/o) G Qt/2,R/2- 



Step 4- It remains to prove the C"-continuity of DW in Q^^2 r/2' 

For two distinct points zi = {xi,yi), Z2 = {x2,y2) £ Qt/2 R/2 '^i^h xi < X2. There 



are two cases: 
Case 1: zi £ R^^- Then 



xi = X2 + ^S, yi = y2 + for some (5, T) G Qi. 



By ()5.33p . we have 



< c. 



Since X2 G (0,?^) and r < 1, the last estimate implies 

\W^{xi,yi) - W^{x2,y2)\ 
((xi - X2)^ + (yi - y2)^)^ 

Case 2: z\ ^ Rz^- Then either |xi— X2I > ^ or \y\—y2\ > Since < X2 < r < 1, 
we find 

\Z1-Z2\ > (y) . 

Thus, using (|5.34p and xi < X2, we obtain 

\W^{zi)-W,{z2)\ ^ \W,{zi)\ + \W^{z2)\ ^^xf + 



\zi — 2:2!" I^l — Z2\°' X2 

Since zi / Z2 are arbitrary points of Q^/2 _r,/2' obtain 

where C depends only on a, r, and the data. The estimate for Wy can be obtained 
similarly. In fact, for these derivatives, we obtain the stronger estimates: For any 

where C depends only on a, r, and the data, but is independent of 5 > and zq. 
Thus, W G C'^'"(Qj',2 D/2) with ||M/^|Li ^,7^+ depending only on the data, 

because r > depends on the data. Moreover, (|5.34p implies that DVF(0, y) = for 
any \y\ < ^. This concludes the proof. □ 
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Theorem 5.6. Let p £ C'^^'^{^l)riC{Q) be the solution of the free boundary problem 

(|2.1ip - (|2.15p in ^4- Then p cannot be across the degenerate sonic boundary 
r 

some- 
Proof. On the contrary, assume that p is across the degenerate sonic boundary 
Tsonic- Then ip = cl — p'^~^ is also across rsonio satisfying the Dirichlet boundary 
problem ()5.3p and (I5.5p ~ (l5.6p in for some small r > 0. Moreover, since ip = ^ 
in the region P1P2OP3, we have Dip = at all (C,?/) £ Tsonic- 

Let (0,7/0) be a point in the relative interior of Fsomc- Then (0, yo) + Q^r C 
if r, R > are sufficiently small. By shifting the coordinates {x,y) — ?• {x,y — yo), 
we can assume (0, j/o) = (0,0) and Q^j^ C Og^. Note that the shifting coordinates 
in the y-direction does not change the expression in (|5.3p . 

Since ip £ C^^^Isq U Lsonic) with Dtp = on Fsonicj reducing r if necessary, we 
obtain 

IV^I < 5x, 

where 5 > is so small that (jS.lip holds in Q^j^- Also, we have ip = c[ — (P' {p) > in 
Q^^. Now we can apply Lemma 5.2 to conclude 

'4>{x,y)> pc2X on Q+is^ 

' 16 

for some p, r > 0. This contradicts the fact that Dip{0,y) = for all y € {—R,R), 
that is, Dip{^,rj) = at any [Cf]) ^ Tsonic- This implies that p cannot be across 
the degenerate sonic boundciry rsonic* 

□ 

We now study more detailed regularity of p near the sonic circle. Prom now on, 
we use a localized version of fig: For given neighborhood A/'(rsomc) of Fsonic and 
e > 0, define 

ne:=nn AA(Fsonic) n{x <e}. 

Since AA(Fsonic) is fixed in the following theorem, we do not specify the dependence 

of on M {T sonic)- 

We first show the regularity part of Theorem 12.11 

Theorem 5.7. Let p be the solution of the free boundary problem (j4.15p - (|4.19p 

established in §^ and satisfy the properties: There exists a neighborhood NiT sonic) 
of T sonic such that, for'ip:=c\- c^{p), 

(a) ip is C^'^ across the part T sonic of the degenerate sonic boundary; 

(b) there exists tDq > so that, in the coordinates (15. ip . 
(5.35) \iP\ < (2ci - ^o)x innn Af{Tsonic)- 
Then we have 

(i) There exists eo > such that ip is C^'" in Q up to Fgonic away from the 
point Pi for any a G (0, 1). That is, for any a G (0, 1) and any (Co;%) G 
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rsonic\-Pi, there exists K < oo depending only on pQ, pi, 7, Eq, a, \\ip\\co,i, 
and d = dist((^0) %)) rsonic) so that 

IIV'lli,„;B^(5o,»;o)nf^,Q/2 - 

(ii) For any (^0,%) £ TsonicVA, 

lim Drip = ci; 

(iii) The limit lim Drib does not exist. 

{5,»?)6n 

The proof is quite similar to the one in [2], which can be achieved by following 
the proof of Theorem 4.2 in [2] step by step. Hence we omit the proof here. 

6. Proof of Theorem 2.1: Global Solutions 

Finally, we show that the solution established above is a global solution indeed, 
valid through the sonic circle rgonici as claimed in Theorem 2.1. 

Since p is only Lipschitz continuous across the sonic circle, we treat the solution 
in the weak sense: For every ^ G C^(r2_), with denoting the region of the left 
state, 

r 2 2 

/ ((c2 - r^)prCr + ^PeCe - -PrC) drde = 0. 

Notice that p is Lipschitz continuous across the sonic circle. Then, due to the Green 
theorem, the integrand is equal to if and only if 

(6.1) [[{c^ - r'^)pr,^pe) ■ 1^] = ^ on Fsonic, 

where the bracket [•] denotes the difference of the quantity between two sides of 
the sonic circle, and v is the normal direction. It is obvious because from the facts 
that {pr,pe) = (— ci,0) up to the sonic circle from the subsonic domain obtained in 
Lemma [531 {pr, pe) = (0, 0) from the supersonic domain and the fact that c^ — = 
on the sonic circle. 

This completes the proof of Theorem 2.1. 

7. Existence and Regularity of Global Solutions of the Nonlinear 

Wave System 

In our main theorem. Theorem 12.11 we have constructed a global solution p of 
the second-order partial differential equation M.lSh in combining this function 
with p = pi in state (1) and p = pQ in state (0). That is, we have obtained the 
global density function p that is piecewise constant in the supersonic region, which 
is Lipschitz continuous across the degenerate sonic boundary Fgonk from 0, to state 

(1)- 
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To recover the momentum components, m and n, we can integrate the second 
and the third equation in (jl.7p , which are equivalent to the equations in (j2.8p in the 
polar coordinates. These can be also written in the radial variable r, 

din 1 Qti 1 

(7-1) — = -p{p)^, -TT- = -P{p)rj, 

or r or r 

and integrated from the boundary of the subsonic region toward the origin. 

Note that we have proved that the limit of Dp does not exist at Pi as r]) in 
fl tends to but \Dc{p)\ has a upper bound. Thus, p{p) is Lipschitz, which 

implies that {m,n) are at least Lipschitz across the sonic circle Tsonic. 

Furthermore, (m, n) have the same regularity as p inside except the origin r = 0. 
However, (m, n) may be multi-valued at the origin r = 0. 

In conclusion, we have 

Theorem 7.1. Let the wedge angle 9^ be between — vr and 0. Then there exists a 
global solution {p, m, n){r, 6) with the free boundary r = r{6), 6 G [6^, Oi], o/ Problem 
2 such that 

(p,m,n) GC2+"(n), peC^m, reC^+''{[9^,9i))nC^^\[9^,9i]), 

and {p,m,n) = {pi,mi,0) in the domain < ^i,r > ri} and (/3o,0,0) in the do- 
main > (,1,1] > r]i}L){r > r{9),9 G [9w,9i]}- Moreover, the solution {p,m,n){r, 9) 
with the free boundary r = r{9) satisfies the following properties: 

(i) p > po on the shock Tshock; that is, the shock Tshock is separated from the 
sonic circle Cq of state (0); 

(ii) The shock Tshock is convex in the self-similar coordinates (C,^)/ 

(iii) The solution {p, m, n) is C^'" up to Psonic CLnd Lipschitz continuous across 
r 

^ some; 

(iv) The Lipschitz regularity of the solution across Psonic (^nd at P\ from the inside 
is optimal; 

(v) The momentum components (m, n) may be multi-valued at the origin. 



Appendix A. Proofs of Lemmas 3.4, 3.8, 3.9, and Proposition 5.3 

In this appendix, we give detailed proofs of Lemmas 13. 4|, [3^81 and l3.91 and Proposi- 
tion [5]3] to make the arguments self-contained. The proofs here are slightly different 
from the ones in [2], [5], and [18], respectively, in order to be adapted for our problem, 
but the main ideas are the similar. 

A.l. Proof of Lemma 13.41 Away from a neighborhood Bd^-^{P2) of P2, the op- 
erator M is oblique. Thus we can apply Theorem 6.30 in [16] to obtain (I3.17P in 
P(d)\{i?rf(j (Pi) U Prfo(P2)}, with a constant C depending on e, ai, VL, do, and Kq. 
Hence we consider only the estimates near P2 in the rest of the proof. 
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For a given solution u to (j3.12P " (j3.13p . we define 

2 

(A.l) V = „ and z = = V (3i(P)DiV. 

l + \\Du\\o ^ 

For do > small enough, O ^ Bdg{P2)- 

Let C be the regularized distance function (from the boundary component Fghock); 
which has the properties: 

1 < ^ < 2, < Co < \DC\ < Cd, \D^C\ < CD(i"^~' 

for two positive constants Co and Cd (see [20]). A smooth approximation to d{P) = 
dist(P, Fshock) is necessary since Fshock has minimal regularity. 

We first construct barrier functions zizg for z on B := -8^0(^2) H fi, by finding a 
suitable positive, increasing function g, g{0) = 0, such that 

1^1 < 9- 

We now derive some estimates for 2 ± / hereafter. 
First, we have 

2 

Diiz + g) = Y^ {PjDrjV + DiPjDjv) + g'Ci, 

that is, 

2 2 
(A.2) Ai^' = Di{z + g)-{Y, D^(3,D,v + g'Q. 

j=i i=i 

We also have 

2 

(A.3) Dij{z + g) = ^{pkDijkV + DjPkDikV + DiPkDjkV + Dij^kDkv) 

k=l 

+9'Cij + g'CiCj- 

In addition, since w satisfies (W2) in Definition 3.2 with a given constant K, we 
obtain the following estimates on the derivatives of aij: 

(A.4) I (a^,) I < WljJ + \alj^u\\Dir,0)w\ < MjJ + \al^,J\\w\\i^' d''^'' < md^'~\ 

(A.5) |L»J^_g)(a^j)l < \(4j,x,x\+'^\(4j,x,u\\D(rfi)W\ + |afj-„^„||L»(r^e)^«P + Kj,nll-C'?r,0)^l 

<l4,x,J+2|af,,.,J|h|ir^ci^-^ + |af,.,,J(||u;||-^M^-i)2 

+ I4,nllkll2""^rf"^"' 
< ?n((i^l-2 + (i27i-2)_ 

Here the subscripts denote the partial derivatives of a^j{r,6) as functions of {x,u). 
The symbol m = m{K) denotes a quantity that depends on the structure of the 
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derivative of a?- and the bound K on w. We absorb the terms which are less 
singular as (i — )■ 0. In a similar way, we obtain the estimates on the derivatives of 

/3 = (/3i,/52). 

Let 72 = min{7i,Q;i}. Then, 

(A.6) \DPi\ < md?^-^, \D'^l3i\ < m{d^^~^ + d^^^'^), 

where m = m{K) > depends on the structure of the derivatives of /3. Here, 
we have used the fact that r', r" , and r'" are bounded by d"^, and d"^"^, 

respectively, as we can apply Lemma 2.8 of [15] to r{9). 

Since l3i{P2,w) = and f32iP2,w) ^ 0, by continuity, we have l32{P2,'w) 7^ 
in Bdo{P2), < do < di, with di depending on and ||r||ci, which is small 

enough. Now we solve the two equations in ()A.2p along with Lv = 0, i.e., 

(A. 7) afjDijV = —{DjaljDiV + bfDiv), 

as a linear system for the three derivatives Dijv. 

The assumption that /32 is bounded away from zero, coupled with the ellipticity 
of L, gives a uniform bound ci(A, A, ||/3||o) on the inverse of the coefficient matrix of 
the linear system. Here we may let A and A be the eigenvalues of afj restricted to 
B, which are bigger than a positive constant depending only on the data and 5 > 0. 
Furthermore, we can use (IA.4|) and (IA.6P to estimate the right-hand sides of (IA.7P 
respectively. We obtain 

(A.8) \D\\ < ci(A, A, \\l3^'\\o){\Diz + /)| + {md^^"^ + W\\o)\Dv\ + fCn). 

This controls the second derivatives of v in terms of |-D(z + /)|. 

Now we proceed to obtain the bounds for z + f. The idea is to find a second- 
order elliptic operator, oi z + g > 0, in B involving the third derivative of v and 
simultaneously to force z + g > on dB, by choice of the function g. We estimate 
these first. Using Lv = 0, \Dv\ < 1, and (jA.Sp . we have 

afjDijkV = - [DkaljDijV + DjafjDikV + bfDikV + Dj^afjOiV + DkbfDiv) 

< (md^2-i + \\b^\\^)\D^v\ + {md''^"^ + md^^^-^ + \\b'\\i)\Dv\ 

< Ci{md^'''~^ + \\b'\\o)\D{z + g)\ + ci{md^^"^ + Wb'Wof 
+Ci{md^^-^ + ||6^||o)5'Cd + ^ ,^^272-2 ^ ||^e||^ 

< C2(^{md^^-^ + l)\D{z + g)\ + {md^^-^ + 1)^ 
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where C2 = C2{A,X,p2, ||/3^^||o, ||6||i,Cd)- Thus, using (|A.7p and makmg the esti- 
mates indicated, we have 

(A.9) atjD,,iz + g) 

< C2m\o(^imd^^~^ + l)\D{z + g)\+{md'^^-^ + l)^ 

+ (m(i^2-i + 1)^' + mci^2-2 ^ ,^^272-2 ^ 

+2Amd^^~^ci (\D{z + g)\+ md^^-^ + ||6||o + g'Co) 
+Am{d''^~^ + d^''^-^) + Ag'\Cij\ + g"dijCiCj 

< Cs(^{md^^-^ + 1)\D{Z + 5)1 + m(i^2-2 ^ ^^2 ^ „^)^272-2 ^ ^^^^72-1(1 + 

Here C3 > is a constant depending on the same parameters as Ci and C2, and the 
terms that are bounded as d — )• have again been omitted. 
Now we define 

(A.IO) Li(z + g) = aljDijiz + g) - c^{md'<^~^ + l)\D{z + g)\. 

Then we have 

Li{z + g)< C^{rnd^^-^ + (m^ + m)d272-2 ^ „,rf72-i(i + ^/)^ + a^'K.jI + g"a%CiCr 
To obtain this estimate, we have assumed that g" < and estimated that 
(A.ll) 9"aljQCj < 9" min {af^QCj) < /A|i^CP < g"Kl 

where we have used the fact that \D(^\ > (^q. 

Now we specify g{() = 5oC^ for any /U < 72. Thus, we have 

/ = goKf^ - 1)C"~' < 2^-'5o^Km - < 0, 

and 

We choose go large enough and d2 G (0, 1) smah enough to obtain 

Li{z + g)<0 inBdo 

for every do < d2- We define do := min{(ii, ^2}- 

Additionahy, ()3.17p is vahd on dB since it holds near Fshock and away from Pi 
and P2. We can choose fo larger if necessary so that (z + g) > on dB. Therefore, 
by the maximum principle, z + g > in B. Thus, z > —g in B. 

Similarly, for z — g{C), it is easy to check that 

aljDij{z-g) 

> -Csl^imd^^-^ + l)\D{z -g)\+ md^^'^ + (m^ + m)d^^^-^ + md^^-^{l + 5')) 
-Ag'\Cij\-g"aljCiCr 
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We define 

(A.12) L2{z - g) := aljD,j{z - g) + c^imd^^-^ + l)\D{z - g)\. 

By a similar argument, we have L2'i/'2 > in B. Therefore, we have 

z < g in B. 
Let F := Li{z + g) < Cd^'^"^. It follows from ([AlO]) that 

aljDij{z + g) = C3{md^^~' + l)\D{z + g)\ + F 

< C3{md^^-^ + l)\D{z + g)\+Cdt'-^. 
Moreover, we can use the Schauder estimates (cf. Lemma 6.20 in [TB]) to obtain 

< C{md^^~^ + l)\D{z + g)\+Cd''-^. 
Using the interpolation inequality, we have 

\\z + gt/ls<Cim). 

Finally, this leads to 

(A.14) \D{z + g)\ < ||(z + 5)||Si_;^';)dA-i < c(m)#-i for d < do. 

We now use (|A.14p in (jA.Sp and drop the lower-order terms to obtain 

\D\\ < ci(|D^i| + md^2-i + g'^ < cd^'-K 
Using Lemma 2.1 in [15], we have 

\\Dv\\^ = \\Dv\\~^^<C{f,)\\Dv\\^^, 

which leads to 

lbl|i+M < c. 

Finally, using the definition of v in (jA.ip . we apply the interpolation inequality with 
small i9 > to obtain 

(A.15) ||n||i+^ < C(l + \\Du\\o) < C(l + + C^\\u\\o) 

and thus (|3.17p holds. Therefore, we obtain the Holder gradient estimate at Fshock 
for the solution u of (j3.12p . This completes the proof. 



A. 2. Proof of Lemma 13.81 For the notational simplicity, throughout the proof, 
we write p = p^'^ . 

We show our claim by contradiction. More precisely, assume that there exists a 
non-empty set B = {X S Fghock ■ c{p, po) ~ ^ > 0} and a point X £ B such that 

niax (c^(p,po) - r"^) = (c^ - = m > 0. 

B 

It is clear that X ^ Pi, P2. Therefore, if X exists, then X £ Fshock \ {-Pij ^2}- Then 
X can be either a local maximum point or a saddle point in U Fghock- We show 
that both cases can not occur, which implies that such X does not exist. 
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First, if X is a local maximum point, then the tangential derivative becomes 

r'{{(?ypr - 2r) + {<?)pe = 0, 
while the normal derivative at X becomes 

(c2)Vr - 2r - r'(c2)Ve > 0, 

which leads to 

(A.16) ((c2)Vr - 2r) (1 + {r'f) > at X 

On the other hand, multiplying c' both sides of Mp = and using the equation due 
to the tangential derivative at X, it follows from (IA.16P that 

2 

1=1 

= /3i(c2)V + /32V'(-(c')V + 2r) 
= (c2)V(/3i-/32r') + 2rr'/32 

= -2(c2)Vrr'(c2 - c^y + 2rr'{3c\r^ - c^) - - r^)) 
< -4rV(c^ - c^) + 6rr'c^(r2 - c^) < 0, 

since we have r'{9) > 0. This is a contradiction. 

Next, if X is a saddle point, then multiplying (c^)' both sides of p = yields 

= (c^YQ^p 

= E-=i«f.(A^(c-2) - (f^iAc-^p) + (#y + M^)i)+nc%, 

where a = ^ and = 6^ — 2r. Thus we can write 

2 
i=l 

where ai = -af^ + ^ and a2 = -a^22 + 7^^- 

Since X is a saddle point, we can construct a barrier function ip so that X = 
{fx, (^x) is a maximum point along the normal direction. 

We define d := Vx — r + r'{6x){9 — Ox) and a set 

W := {(r, e)en:d>0}n {(r, 9) e Q : c'^ - > m}. 

Note that, since X is a saddle point, there are interior points of the set W, especially 
in the inward normal direction. We note that, since c^(p,po) > r'^ + m in W and 
c{p) > c{p,po) for p > pq, there exists a constant eo > such that 

c^(/o) - > eo > in W. 
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Set u = — r'^ — m. Then 

2 

= ^(af,Ai(c2) + a.|A(c2)P) +Kc')r 
i=l 

2 

= [aliDii{u + + m) + ai\Di{u + + m)p) + b{u + + m)^ 

7=1 
2 

= ^ {afiDiiU + Oil A/up) + (4rai + + 2ali + 4r2ai + 2r6^, 
1=1 

where 6o 4air + 6^ and /o := 2af + Ar'^ai + 2r6'^. We now let 

:= ^(ew« _ 1), ;Uo>0. 

Then we have 

2 

= {atiDiiU + ai\Diu\'^) + bpUr + /o 

2 

Hence, we have 

2 2 

5^ afiA*^« + bowr + e^o^/o > 6"^°" J^(af,Mo - a.)| A.^^l' > 0, 
1=1 1=1 

by choosing /^o = max{aj}/eo, where af^,a22 > cq > in W, and /^o and eo are 
independent of e. 

Thus, we find il^{d) with -;/''> and ip" < satisfying 

2 

Li^ = a,, A*w^ + boWr + e'^°"/o < cq^" + + /i < 0, 

where 6i = niax(4air + b), fi = maxe'^°"(2a^° + 4r^ai + 2r6^)+, and e < eo- The 

w w 
solution to the equation: 

eo^" + biij' + A = 

is 

6i 1 - e-^i'^i/^o 6i ' 
which satisfies the boundary condition: 

V^(O) = 0, tP{di) = mo, 

with mo = where Umax = = niax(c^ — — m) and di > 0. 
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Hence, in the set W, we have 

Liiw-i)) > 0. 

Thus, by the weak maximum principle, we obtain 

sup(t(; — Tp) = sup(tt; — ip)'^ < = {w — i/j){X). 
w aw 

Thus, X is the local maximum point of u = w—ipmWso that the normal derivatives 
at X becomes 

(A.17) y^-r'{d)vg = —{X)>0, 

while the tangential derivatives at X along Tshock becomes 

dv 

(A.18) r'vr + vg = r'{wr - ip'dr) + {we - i^'dg) = —{X)= 0. 

Hence, we obtain 

Vr{l + {r'f) >0, 

that is, 

(A. 19) Wr - Ip'dr = Vr{X) > 0. 

Now, since p satisfies Mp = on Tshock with given r'{9) > 0, multiplying e^o^c' 
throughout Mp = 0, it follows from ()A.18p that 

2 
i=l 

= /3iK + 2e'^°V)+/32i/;9 

= 2e^o"r/3i + (/3i - /32r')u;, + f32{i^'de + /SaV'V'dr). 
At X, noting that n = and using (lA.lSp . we obtain 

> 2r/3i - iP'pdr{X) + hiip'de + M'r'dr){X) 

= 2r/3i+V'(0)(/3i-r'/32) 

= 2r/3i +V'(0)Ai. 
Thus, we show that, if ip satisfies 

'^^ ' ~ eo(l - e-^i'^i/^o) bi p 

so that the last inequality becomes negative, which is a contradiction, where we have 
used the fact that ^ < 0. 

The second order Taylor series expansion of ^-^^^,1-/^0 

61 . h 



eo 2eg 
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for some < (i2 < di- Hence we have 

where = 0{di) = Ae-^^'^^/^°di. 

Now, at X, recah (jA.lSp and notice that 

r'wr + wq = r'tp'dr + ip'dg = 0, 

so that 

(A.20) r'(c^)Vr - 2rr' = r'wr = -we = -{c'^)'pe- 

Also, from Mp = 0, we have 

r{c)pr = -r{c) —pe = (-1 + -^)(c ) pe- 
Pi Pi 

Thus, from (|A.20p . we have 



2w _ 2rV/3i ^ _ 2r/3i 



(A.21) we = {c'Ype = w, 

P P 

Notice that, for small di > 0, we can write rriQ = w{P) at some point in W and 
mQ = w{P) < w{X) + Dw{P)-{P-X) + \\w\\i+c,,wd\^'' 

where mi = \\w\\i^a,w- Thus, for all di, we obtain 



where we recall = 0{di). Then we have 



2r/3i + V'(0)/i < 2r/3i + + midf + {^^ 

I — v\ p bi 

Therefore, we have 

2r/3i + ^ = ^r'(0)(c2(r2 - c^) - 3c\c^ - r^)) < for t9 < 1, 

I — V 1 — V 

where we have used the fact r{X) < c{X). Thus, for sufficiently small di and 
1? = 0{di), we have 

2rPi+ij'{0)p < 0, 

which is a contradiction. 

Therefore, there is no such X, which implies that the set 5 = 0. This completes 
the proof. 
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A. 3. Proof of Lemma 13.91 For the notational simplicity, throughout the proof, 
we write p = p^'^ . To prove the monotonicity, we argue by contradiction. 

Let us first examine the C^-function p restricted to Lshock- Since r'{6) > 0, p is 
now a function of a single variable, i.e., the second component of a point P = {r{0), 6) 
on Lshock) without confusion. We can order the points along Fghock by 6 and refer 
to the intervals along Lshock by the label. Then the lack of monotonicity implies 
that there exist points 0i and G2 on Lshocki with P2 < ©i < ©2 < ^i, at which 
p(0i) > /o(02). Then we immediately deduce that 

1. In (^2)02)) there exists C with p{C) = max p; 

[-P2,©2] 

2. In (C, Pi), there exists D with p{D) = min p. 

We want to identify the points C and D on Pshock with C < D such that 

(i) piP2)<P<p{C) on [P2,C]; 

(ii) p{C)>p>piD) on [C,D]; 

(iii) p{D)<p<p{Pi) on [D,Pi]. 

Now, property (ii) may not hold with C = C because p{C) is the maximum value of 
p only at the interval [P2) D], and we may have -D > ©2. Then, if there is a point in 
{P2-, ©2) at which p > p(C'), we let C to be the point. Otherwise, we choose C = C. 
Thus, all the three properties hold. 

Now we look at the function p in Vt. The idea is to partition into three subdo- 
mains by two curves Tq and Fd from C and D to the points A and B respectively 
on Fq, in such a way p{A) > p{B) that we can deduce that there is a point m on 
Fq at which p obtains a maximum on either the subdomain Q.a or the domain fi^, 
thus violating the Hopf maximum principle. This is also the case even if it happens 
to be the origin O. It suffices to show that p{rn) is the maximum value of p on the 
boundary of Q,a ov VLb- 

We now construct the Lipschitz curves on which p has certain monotone property. 
That is, 

(A.22) p{A)>p>p{C)-iionTc, p{A)>p{C), 

and 

(A.23) p{B) < p < p{D) + p ouFd, p{B)<p{D), 

for certain number p > 0. We specify 

p = \ min{p(C) - p(D), pi - p{C),p{D) - p}. 

Since p G C°(r2), we have 

\p{Xi) - p{X2)\ < M\Xi-X2r 
for some M > and Xi,X2 € ^2. Now, on any ball with radius r > 0, 

Osc(/3) < 2Mr°. 
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Let R = (sir)""- We have 

Oscs^nn(p) < 

Now we construct Tc as follows (as in Fig [5]): In Bii{C) Pi 0, piC) can not be the 
maximum value of p. Thus, there exists a point in dBii{C) R where p > p{C). 
Let Xi be a point at which p attains its maximum value in Bji{C). Since the first 
segment of Tc is a straight line from C to Xi , we have 

piXi) > p{C), 

and, on the segment, 

piX)>piC)-p, p{X)<p{Xi). 




Figure 5. Hypothetical Curves 



Now we continue inductively, forming a sequence of the line segments with corners 
at {Xi} (take Xq = C), along which p{X) > p{C) - p and p{Xi) < p{X2) < ■ ■ ■ . 
To show this fact, we let 

f), = ne\{jl^BR{Xi)}. 

Then Xj G 9%. 

Consider Bii[Xi). Note that p{Xj) is the largest value of p on the part of Bpi{Xi) 
inside the complement of VLj. However, p{Xj) is less than the maximum value of p on 
Bpi{Xj), by the mean value property. Hence, there is a point Xj+i G dBii{Xi) n Q,j 
at which p attains its maximum value in Bji{Xi). Again, along the straight line 
from Xj to Xj^i, we have 

p > p{X,) -p> p{C) - p. 



Now, 



dist(Xi_i,fii) = R 
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and 

Qj c c • • • c ill, 

so that, ior k > i + 1, 

dist(Xi,Ofe) > R. 

Since G 5^7^, the estimate 

dist(Xj+i, Xj) >R for i < j 

follows. Hence, 

dist{X j,Xi) > R for i ^ j. 

Since the domain Q is finite, this process must end at finite steps when we reach 
a point Xl = B £ d^l. By construction, Tc satisfies the properties in ()A.22p . 
Similarly, we construct Td, with termination point A S dQ. 

We now locate A and B. We note that the two curves cannot cross each other. 
This is because, at every point on Tc, 

p > p{C) - fi > piD) + fi; 

while, at every point on Fd, we have 

p<p{D) + fi. 

Furthermore, Tq cannot terminate at Fsonic where p > pi — p > p{D) + p. For the 
same reason, it can not come back to Fs^ock in [-P21C] or [C,D] where p < p{C). 
Finally, A cannot lie in the segment [^2, C] of T shock-, because this would trap F^ in 
a region where p > p{D) (or more simply, this would contradict with the fact that 
D is not a local minimum in Q). Hence, A has to end on Fq. 

Similarly, B cannot lie on Fghock where p > p{D) in the interval [D,Pi) and must 
lie on Fq (see Fig [5]) . 

Now we reach to our final contradiction. Since p{A) is larger than p and p{B), 
there is a point m along the boundary P2OB at which p attains a maximum. Assume 
first that m is not the origin, then m can not be a local maximum for the domain il. 
by the Hopf lemma. However, along the entire boundary of the domain P2CDBP2, 
p < /o("^)) which implies that it is a maximum. This is a contradiction. Now, if 
m coincides with O, the similar minimum point X resembling B can not coincide 
with O. We can find that there is no place for such X either. Thus this is also a 
contradiction. We conclude that p is monotone along Fghock from P2 to Pi. 

A. 4. Proof of Proposition 15. 3L The proof is divided into two steps. 

In the proof below, all the constants depend only on the data, i.e., A'^, ci, f, R, 
j3, and inf il), unless otherwise is stated. 

Q+^n{x>f/2} 

Step 1: We prove that there exist ai E (0, \) and ri > such that, li W £ 
CiOfji) n C^iQt,R) satisfies (I5:20]) ~ (l5:22]l . then 

(A.24) W{x, y) < "1(^-/^1) ^1+^ in , 
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whenever a S (0, ai], r G (0, ri], and /^i < min{/i, i}. 
By Lemma 5.2, 



(A.25) W{x,y)<cx'y\- iix)x in Q 

where rg depends only on the data. 

Fix yo with |7/o| < We now prove that 

yo) < £lii:ii^xi+° for x G (0, r) . 



As in [2], without loss of generality, we may assume that R = 2 and yo = 0. 

§y) for all {x,y) G Qt^. 



Otherwise, we set W{x,y) = W{x,yo + ^y) for all G Q/g- Then G 



C(Qt^) n C2(Qt^) satisfies ([5:20]) - (f5:22]) . with some modified constants A^, and 
Oi, depending only on the corresponding quantities in the original equation and on 
R. We conclude that, without loss of generality, we can assume that yo = and 
R = 2. That is, it suffices to prove that 

(A.26) W{x, 0) < for x G (0, r), 

for some r G (0,ro) and a G (0, ai), under the assumptions that ()5.20p - ()5.22|) hold 
in Qt^ and (1X25]) holds in Q+ 2- 
For any given r G (0,ro), let 

(A.27) Air = ci{l-ni), Si = ci(l-^i), 

(A.28) v = Aix^+''{l-y^) + Bixy^. 

Since (fOT]) holds on Qt^nlx = 0} and (|X25]) holds in 2' then, for all x G (0, r) 
and \y\ < 1, we obtain 

(v{0,y)=0 = W{0,y), 

vir, y) = {Air^{l - y^) + Biy^)r = ci(l - /ii) > Ty(r, y), 

v{x, ±1) = ci(l - > W{x, ±1). 

Thus, we have 

(A.29) W <v on SQ+^p 

From ()5.20p . we obtain 

„2/ 



In order to rewrite the right-hand side in a convenient form, we write the term Vyy 

Jyy -r y Uyy . 

\ _ „ „an ,.2\ T I „,2 



in the expression of C2V as (1 — y'^)vyy + y'^Vyy. Then a direct calculation yields 



C2V - (C1O2 + cfOs) = cix"(l - y')Ji + y'J2, 
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where 

= Q,(a + l)ci^i + a(a + l)^?x"(l - y^) + a(a + - (a + llci^i 

X 

+{a + 1)(02 + 2C1O3) + (a + ifAlx^ - [a + ifAlx^O'^x'^ 

+2(a + l)AiSix"03y^ - 2^ix + 2^12;^-" + 04(-2Aix + 2Bix^-") 

-C2— - 2ci — , 

and 

J2 = a(a + l)AiSix"-2ciBi + BiV + 2(a + l)^iBix"(l-y2)_52Q3y2 

+04(-2ylix°+i + 2Bix) + 25ix"+i + 4f- - O5) + Bixf. 

V(7-l)cf / 

By (fS^Sj) and (|X27l) . we obtain 

(A.30) Ji < (a - l)(a + l)ciAi + 2C(2 - //)r^-", 

(A.31) J2 < -Bici + 5? + Cr°. 

Moreover, 
(A.32) 

C2{v -W)- v^v - W) 

= (cix + W + Oi){v - - (ci - O2 - 2ci03)(^ - 

+(1 - 0^)(y + VF),(7; - + (1 + 0^\v - W)yy 
-{j:i^^-0,){v + WUv-W). 

< cix°((a - l)(a + l)ciAi + 2C(2 - ^f)ri-")(l - y^) + (-^ici + 5^ + Cr")?/^ 
-(1 + a)a^i(^ixi+"(l - y^) + Bixy^){l - y^) 

+ (1 + a)aci{l - fii)Aix''{l - y^) 

< cix°((2a - l)(a + l)ciAi + 2C(2 - /i)ri-") (1 - y^) 

+ (-SlCi+B2 + C7r")y2 

= cix"(l-y2)/^ + y2j2, 
where we have used the fact that 

v^AW -v) = a{a + l)Aix"-i(l - y^){W - v) < {1 + a)aci{l - - y^). 

Then we can choose ai > 0, depending only on ^1, so that, if < a < ai, 

(A.33) (2a - l)(a + l)ciyli < -y. 

Such a choice of ai > is possible because we have the strict inequality in ()A.33P 
when Q = 0, and the left-hand side is an increasing function of a > 0. Now, choosing 
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ri > SO that 

(A.34) n<min{(i;i)i(^l^)iro} 

is satisfied, we use ()A.30P ~ ()A.34p to obtain 

(A.35) Ii < 0, /2 < in Q+i- 

Then, by ()A.32p . we obtain 

(A.36) C2{v-W)-v^^{v-W)<{) inQ+ 

whenever r G (0,'ri] and a G (0,ai]. Then v — W satisfies the maximum principle 
in Q^i, which imphes 

(A.37) W <v, in Q+i. 

In particular, using ()A.27p and ()A.28P with y = 0, we arrive at ()A.24p . 

Step 2. We now establish Proposition 5.3. As argued before, without loss of 
generality, we may assume that R = 2 and it suffices to show that 

(A.38) W{x,0) < ^a;i+" for x £ [0,r]. 

By Step 1, it suffices to prove that (|5.23p holds for the case a > ai. Fix any 
a G (cti, 1) and set the following comparison function: 

ri r 

By Step 1, 

(A.40) W<v on ag+i for r G (0,ri]. 

As in the proof of Step 1, we write 

r.,iv -W) = {1 + a fil -J^'K -il - yVi + (1 + ai)^i^^^x-y2j2, 



r-^ r 



where 



- Cl(l — Ul) 2n ci(1 — Ul) 2 Oi 

oi J.^^ J. 



+0, + 2c,0, + (1 + a)(l + Osfj^^x'^il - y^) 

+2(1 + 03)(1 + ai)^^^^x°iy2 _ ^ + ^ i^:^ x 

r" a + 1 a + 1 

a + 11 — ^1 X x 
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and 



r" "1 x 



1 rfxi-^i O2 , O3. 



1 + ai 1 — ^1 2; X 

It is easy to see that Ji, J2 < 0. By continuity, we can choose r > sufficiently 
smah, depending only on N, ci, f, ii, and ai. Thus, we obtain 

C,{v-W)-v,,{v-W) < {l + afj^-^x^{l-y'){{a-l)c, + Cr^^) 



+(1 + a, f'^\^'K "^y'{{a, - l)c, + Cr"^) 



1 

where we have used the result in Lemma 5.3 that W < It is easy to check 

that 

/:2V - C2W - v^^{v -W) <0. 
As the proof in Lemma 5.3, it is easy to prove that (j5.23p holds with 

ci(l - m) 



A 

This complete the proof. 
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